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Abstract 

A description of dislocations and disclinations defects in terms of Riemann- 
Cartan geometry is given, with the curvature and torsion tensors being interpreted 
as the surface densities of the Frank and Burgers vectors, respectively. A new 
free energy expression describing the static distribution of defects is presented, and 
equations of nonlinear elasticity theory are used to specify the coordinate system. 
Application of the Lorentz gauge leads to equations for the principal chiral SO(3)- 
field. In the defect-free case, the geometric model reduces to elasticity theory for 
the displacement vector field and to a principal chiral SO(3)-field model for the spin 
structure. As illustrated by the example of a wedge dislocation, elasticity theory 
reproduces only the linear approximation of the geometric theory of defects. It is 
shown that the equations of asymmetric elasticity theory for the Cosserat media can 
also be naturally incorporated into the geometric theory as the gauge conditions. As 
an application of the theory, phonon scattering on a wedge dislocation is considered. 
The energy spectrum of impurity in the field of a wedge dislocation is also discussed. 

1 Introduction 

Many solids have a crystalline structure. However, ideal crystals are absent in nature, 
and most of their physical properties, such as plasticity, melting, growth, etc., are defined 
by defects of the crystalline structure. Therefore, a study of defects is a topical scientific 
question of importance for applications in the first place. A broad experimental and the- 
oretical investigations of defects in crystals started in the 1930s and continues nowadays. 
At present, a fundamental theory of defects is absent in spite of the existence of dozens 
of monographs and thousands of articles. 

One of the most promising approaches to the theory of defects is based on Riemann- 
Cartan geometry, which involves nontrivial metric and torsion. In this approach, a crystal 
is considered as a continuous elastic medium with a spin structure. If the displacement 
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vector field is a smooth function, tlien tliere are only elastic stresses corresponding to 
diffeomorphisms of the Euclidean space. If the displacement vector field has disconti- 
nuities, then we are saying that there are defects in the elastic structure. Defects in the 
elastic structure are called dislocations and lead to the appearance of nontrivial geometry. 
Precisely, they correspond to a nonzero torsion tensor, equal to the surface density of the 
Burgers vector. 

The idea to relate torsion to dislocations appeared in the 1950s [1-4]. This approach 
is still being successfully developed (note reviews [5-11]), and is often called the gauge 
theory of dislocations. A similar approach is also being developed in gravity ^2j. It is 
interesting to note that E Cartan introduced torsion in geometry [T^ having the analogy 
with mechanics of elastic media in mind. 

The gauge approach to the theory of defects is being developed successfully, and 
interesting results are being obtained in this way [14-17]. We note in this connection 
two respects in which the approach proposed below is essentially different. In the gauge 
models of dislocations based on the translational group or on the semidirect product of the 
rotational group with translations, one usually chooses the distortion and displacement 
fields as independent variables. It is always possible to fix the invariance under local 
translations such that the displacement field becomes zero because it transforms by simple 
translation under the action of the translational group. In this sense, the displacement 
field is the gauge parameter of local translations, and physical observables are independent 
of it in gauge-invariant models. 

The other disadvantage of the gauge approach is the equations of equilibrium. Ein- 
stein type equations are usually considered for distortion or vielbein, with the right hand 
side depending on the stress tensor. This appears unacceptable from the physical point of 
view because of the following reason. Consider, for example, one straight edge dislocation. 
In this case, the elastic stress field differs from zero everywhere. Then the torsion tensor 
(or curvature) is also nontrivial in the whole space due to the equations of equilibrium. 
This is wrong from our point of view. Indeed, we can consider an arbitrary domain of 
medium outside the cutting surface and look at the creation process for an edge disloca- 
tion. The chosen domain was a part of the Euclidean space with identically zero torsion 
and curvature before the defect creation. It is clear that torsion and curvature remain 
zero because the process of dislocation formation is a diffeomorphism for the considered 
domain. In addition, the cutting surface may be chosen arbitrary for the defect creation, 
leaving the dislocation axis unchanged. Then it follows that torsion and curvature must 
be zero everywhere except at the axis of dislocation. In other words, the elasticity stress 
tensor can not be the source of dislocations. To avoid the apparent contradiction, we pro- 
pose a radical way out: we do not use the displacement field as an independent variable at 
all. This does not mean that the displacement field does not exist in real crystals. In the 
proposed approach, the displacement field exists and can be computed in those regions of 
medium that do not contain cores of dislocations. In this case, it satisfies the equations 
of nonlinear elasticity theory. 

The proposed geometric approach allows considering other defects that do not relate 
directly to defects of elastic media. 

The intensive investigations of other defects were conducted in parallel with the study 
of dislocations. The point is that many solids have not only elastic properties but also a 
spin structure. For example, there are ferromagnets, liquid crystals, spin glasses. In this 
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case, there are defects in the spin structure which are called disclinations JHl- They arise 
when the director field has discontinuities. The presence of disclinations is also connected 
to nontrivial geometry. Namely, the curvature tensor equals the surface density of the 
Frank vector. The gauge approach based on the rotational group §0(3) was also used 
for describing disclinations [12]. SO(3)-gauge models of spin glasses with defects were 
considered in [201 121] • 

The geometric theory of static distribution of defects which describes both types of 
defects - dislocation and disclinations - from a single standpoint was proposed in [22J. In 
contrast to other approaches, it involves the vielbein and SO(3) connection as the only 
independent variables. The torsion and curvature tensors have direct physical meaning as 
the surface densities of dislocations and disclinations, respectively. Covariant equations 
of equilibrium for the vielbein and SO(3) connection similar to those in a gravity model 
with torsion are postulated. To define the solution uniquely, we must fix the coordinate 
system (fix the gauge) because any solution of the equations of equilibrium is defined up 
to general coordinate transformations and local §0(3) rotations. The elastic gauge for 
the vielbein fl^ and Lorentz gauge for the SO(3) connection 2? were proposed recently. 
We stress that the notions of a displacement vector and rotational angle are completely 
absent in our approach. These notions can be introduced only in those domains where 
defects are absent. In this case, equations for vielbein and §0(3) connection are identically 
satisfied, the elastic gauge reduces to the equations of nonlinear elasticity theory for the 
displacement vector, and the Lorentz gauge leads to the equations for the principal chiral 
SO(3) field. In other words, to fix the coordinate system, we choose two fundamental 
models: the elasticity theory and the principal chiral field model. 

To show the advantages of the geometric approach and to compare it with the elasticity 
theory, we consider in detail a wedge dislocation in the frameworks of the elasticity theory 
and the proposed geometric model. We show that the explicit expression for the metric in 
the geometric approach is simpler and coincides with the induced metric obtained within 
the elasticity theory only for small relative deformations. 

As an application of the geometric theory of defects, we consider two examples in the 
last sections of the present review. First, we solve the problem of phonon scattering on 
a wedge dislocation. The problem of phonon scattering is reduced to the integration of 
equations for extremals for the metric describing a wedge dislocation (because phonons 
move along extremals in the eikonal approximation). As a second application, we consider 
the quantum mechanical problem of impurity or vacancy motion inside a cylinder whose 
axis coincides with a wedge dislocation. The wave functions and energetic spectrum of 
the impurity are found explicitly. 

The presence of defects results in nontrivial Riemann-Cartan geometry. This means 
that for describing the phenomena that relate ingeniously to elastic media, we must make 
changes in the corresponding equations. For example, if a phonon propagation in an 
ideal crystal is described by the wave equation, then the presence of defects is easily 
taken into account. For this, the fiat Euclidean metric has to be replaced by a nontrivial 
metric describing the distribution of defects. The same substitution must be made in the 
Schrodinger equation to describe other quantum effects. It is shown nowadays that the 
presence of defects essentially infiuences physical phenomena. The Schrodinger equation 
in the presence of dislocations was considered in [25-46] for different problems. Problems 
related to the wave or Laplace equations were considered in [47-54]. The influence of the 
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nontrivial metric related to the presence of defects was investigated in electrodynamics 
and hydrodynamics jSnj- Scattering of phonons on straight parallel dislocations was 
studied in [57-59]. 

Another approach to the theory of defects based on affine geometry with nonzero 
nonmetricity tensor was considered recently in pUj . 

2 Elastic deformations 

We consider infinite three dimensional elastic media. We suppose that the undeformed 
medium in the defect-free case is invariant under translations and rotations in some coor- 
dinate system. Then, the medium in this coordinate system y*, i = 1,2, 3, is described by 
the Euclidean metric 6ij = diag (+ ++), and the system of coordinates is called Cartesian. 
Thus, in the undeformed state, we have the Euclidean space with a given Cartesian 
coordinate system. We also assume that torsion (see the Appendix) in the medium equals 
zero. 

Let a point of the medium has coordinates in the ground state. After deformation, 
this point has the coordinates 



dinates x that are coordinates of points of the medium after the deformation and cover 
the whole Euclidean space M.^. In the presence of dislocations, the coordinates do not 
cover the whole M.^ in the general case because part of the media may be removed or, 
conversely, added. Therefore, the system of coordinates related to points of the medium 
after an elastic deformation and defect creation is more preferable. 

In the linear elasticity theory, relative deformations are assumed to be small {djU^ -C 
1). The functions = u^{y{x)) are then components of a vector field that is called 

the displacement vector field and is the basic variable in elasticity theory. 

In the absence of defects, we assume that the displacement field is a smooth vector 
field in the Euclidean space R^. The presence of discontinuities and singularities of the 
displacement field is interpreted as a presence of defects in elastic media. 

In what follows, we consider only static deformations with the displacement field 
independent of time. Then the basic equations of equilibrium for small deformations are 




(1) 



Figure 1: Elastic deformations 




in the initial coordinate system, see Fig. ^ 
The inverse notation is used in the elasticity 
theory. One usually writes — > = x^+u\x). 
These are equivalent because both coordinate 
systems and cover the whole R^. However, 
in the theory of defects considered in the next 
sections, the situation is different. Generally, 
the elastic medium fills the whole Euclidean 
space only in the final state. Here and in what 
follows, we assume that fields depend on coor- 
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(see, e.g., jHH) 



d,<y^' + f = 0, (2) 
a'^ = XS'^ek" + 2fxe'^, (3) 

where a^^ is the stress tensor, which is assumed to be symmetric. The tensor of small 
deformations Sij is given by the symmetrized partial derivative of the displacement vector: 

e^j = ^{diUj + djUi). (4) 

Lowering and raising of the Latin indices is performed with the Euclidean metric 6ij and 
its inverse 6^^ . The letters A and /i denote constants characterizing elastic properties of 
media and are called Lame coefficients. Functions describe the total density of 

nonelastic forces inside the medium. We assume in what follows that such forces are 
absent: /-'(x) = 0. Equation Q is Newton's law, and Eqn (jS)) is Hook's law relating 
stresses to deformations. 

In a Cartesian coordinate system and for small deformations, the difference between 
upper and lower indices disappears because raising and lowering of indices is performed 
with the help of the Euclidean metric. One usually forgets about this difference due 
to this reason, and this is fully justified. But in the presence of defects, the notion of 
Cartesian coordinate system and Euclidean metric is absent, and the indices are raised 
and lowered with the help of Riemannian metric. Therefore, we distinguish the upper 
and lower indices as is accepted in differential geometry, having the transition to elastic 
media with defects in mind. 

The main problem in the linear elasticity theory is the solution of the second-order 
equations for the displacement vector that arise after substitution of Q into with some 
boundary conditions. Many known solutions are in good agreement with experiment. 
Therefore, one may say that equations (jH)), have a solid experimental background. 

We now look at the elastic deformations from the standpoint of differential geometry. 
From the mathematical standpoint, map is itself a diffeomorphism of the Euclidean 
space M^. The Euclidean metric Sij is then induced by the map — > x*. It means that 
in the deformed state, the metric in the linear approximation is given by 

9ij{x) = ^-[^^ki ~ Sij - diUj - djUi = 5ij - 2eij, (5) 

i.e., is defined by the tensor of small deformations (jH). We note that in the linear approx- 
imation, eij{x) = eijiy) and duj/dx'' = duj/dy\ 

In Riemannian geometry, the metric uniquely defines the Levi-Civita connection 
rjj'^(x) (Christoffel's symbols), Eqn pi5|) . We can compute curvature tensor ()119|) for 
these symbols. This tensor equals identically zero, Rijk{x) = 0, because the curvature 
of the Euclidean space is zero, and the map — > is a diffeomorphism. The torsion 
tensor is equal to zero for the same reason. Thus, an elastic deformation of the medium 
corresponds to the trivial Riemann-Cartan geometry, with zero curvature and torsion 
tensors. 

The physical interpretation of metric © is as follows. The external observer fixes 
Cartesian coordinate system corresponding to the ground undeformed state of the medium. 
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The medium is then deformed, and external observer discovers that the metric becomes 
nontrivial in this coordinate system. If we assume that elastic perturbations in the medium 
(phonons) propagate along extremals (lines of minimal length), then their trajectories in 
the deformed medium are defined by Eqns (|118|) . Trajectories of phonons are now not 
straight lines because the Christoffel's symbols are nontrivial {Tj^^ 7^ 0). In this sense, 
metric (0) is observable. Here, we see the essential role of the Cartesian coordinate system 
defined by the undeformed state, with which the measurement process is connected. 

We assume that the metric gij{x) given in the Cartesian coordinates corresponds to 
some state of elastic media without defects. The displacement vector is then defined 
by the system of equations and its integrability conditions are the equality of the 
curvature tensor to zero, in accordance with theorem El in the Appendix. In the lin- 
ear approximation, these conditions are known in elasticity theory as the Saint-Venant 
integrability conditions. 

We make a remark that is important for the following consideration. For appropriate 
boundary conditions, the solution of the elasticity theory equations © is unique. 
From the geometric standpoint, this means that elasticity theory fixes diffeomorphisms. 
This fact is used in the geometric theory of defects. Equations of nonlinear elasticity 
theory written in terms of the metric or vielbein are used for fixing the coordinate system. 

3 Dislocations 

We start with the description of linear dislocations in elastic media (see, e.g., jHH E^)- 
The simplest and most undespread examples of linear dislocations are shown in Fig. |21 We 
cut the medium along the half-plane x'^ = 0, > 0, move the upper part of the medium 
located over the cut > 0, > by the vector b towards the dislocation axis x^, 
and glue the cutting surfaces. The vector b is called the Burgers vector. In the general 
case, the Burgers vector may not be constant on the cut. For the edge dislocation, it 
varies from zero to some constant value b as it moves from the dislocation axis. After the 
gluing, the media comes to the equilibrium state called the edge dislocation, see Fig. 
If the Burgers vector is parallel to the dislocation line, it is called the screw dislocation 
(Fig.H^). 

The same dislocation can be made in different ways. For example, if the Burgers 
vector is perpendicular to the cutting plane and directed from it in the considered cases, 
then the produced cavity must be filled with medium before gluing. It is easy to imagine 
that the edge dislocation is also obtained as a result, but rotated by the angle 7i/2 around 
the x^ axis. This example shows that a dislocation is characterized not by the cutting 
surface but by the dislocation line and the Burgers vector. 

From the topological standpoint, the medium containing several dislocations or even 
the infinite number of them is the Euclidean space M^. In contrast to the case of elastic 
deformations, the displacement vector in the presence of dislocations is no longer a smooth 
function because of the presence of cutting surfaces. At the same time, we assume that 
partial derivatives of the displacement vector dju^ (the distortion tensor) are smooth 
functions on the cutting surface. This assumption is justified physically because these 
derivatives define deformation tensor Q. In its turn, partial derivatives of the deformation 
tensor must exist and be smooth functions in the equilibrium state everywhere except, 
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Figure 2: Straight linear dislocations, (a) The edge dislocation. The Burgers vector b is 
perpendicular to the dislocation line, {b) The screw dislocation. The Burgers vector b is 
parallel to the dislocation line. 

possibly, the dislocation axis, because otherwise equations of equilibrium Q have no 
meaning. We assume that the metric and vielbein are smooth functions everywhere in 
except, may be, dislocation axes, because the deformation tensor defines the induced 
metric (0). 

The main idea of the geometric approach amounts to the following. To describe 
single dislocations in the framework of elasticity theory we must solve equations for the 
displacement vector with some boundary conditions on the cuts. This is possible for small 
number of dislocations. But, with an increasing number of dislocations, the boundary 
conditions become so complicated that the solution of the problem becomes unrealistic. 
Besides, one and the same dislocation can be created by different cuts which leads to 
an ambiguity in the displacement vector field. Another shortcoming of this approach is 
that it cannot be applied to the description of a continuous distribution of dislocations 
because the displacement vector field does not exist in this case at all because it must have 
discontinuities at every point. In the geometric approach, the basic variable is the vielbein 
which by assumption is a smooth function everywhere except, possibly, dislocation axes. 
We postulate new equations for the vielbein (see sectional). In the geometric approach, the 
transition from a finite number of dislocations to their continuous distribution is simple 
and natural. In that way, the smoothing of singularities occurs on dislocation axes in 
analogy with smoothing of mass distribution for point particles in passing to continuous 
media. 

We now develop the formalism of the geometric approach. In a general defect-present 
case, we do not have a preferred Cartesian coordinate frame in the equilibrium because 
there is no symmetry. Therefore, we consider arbitrary coordinates x^, /i = 1, 2, 3, in R^. 
We use Greek letters for coordinates allowing arbitrary coordinate changes. Then the 
Burgers vector can be expressed as the integral of the displacement vector 





where C is a closed contour surrounding the dislocation axis. Fig. |31 
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This integral is invariant under arbitrary coordinate transformations 



(x) 



and covariant under global SO(3)-rotations of y^. Here, components of the displacement 
vector field are considered with respect to the orthonormal basis in the tangent 

space, u = u'Ci. If components of the displacement vector field are considered with 
respect to the coordinate basis u = u^d^, the invariance of the integral (jH)) under general 
coordinate changes is violated. 

In the geometric approach, we introduce 
new independent variable - the vielbein - in- 
stead of partial derivatives 



[X 




outside the cut, 
on the cut. 



(7) 



The vielbein is a smooth function on the cut by 
construction. We note that if the vielbein was 
simply defined as the partial derivative d^y^, 
then it would have the 5-function singularity 
on the cut because functions y\x) have a jump. 
The Burgers vector can be expressed through 
the integral over a surface S having contour C 
as the boundary 




dx^e,, 



c 



dx^'Adx'^id, 



d, 



Figure 3: Section of the media with the 
edge dislocation. C is the integration 
contour for the Burgers vector b. 



where dx^ /\dx^ is the surface element. As a consequence of the definition of the vielbein in 
((7j), the integrand is equal to zero everywhere except at the dislocation axis. For the edge 
dislocation with constant Burgers vector, the integrand has a 5-function singularity at the 
origin. The criterion for the presence of a dislocation is a violation of the integrability 
conditions for the system of equations d^y^ = e 



d e * 



- d,e; ^ 0. 



(9) 



we 



If dislocations are absent, then the functions y^{x) exist and define transformation to a 
Cartesian coordinates frame. 

In the geometric theory of defects, the field e^* is identified with the vielbein. Next 
compare the integrand in (jH)) with the expression for the torsion in Cartan variables 
They differ only by terms containing the §0(3) connection. This is the ground for the 
introduction of the following postulate. In the geometric theory of defects, the Burgers 
vector corresponding to a surface S is defined by the integral of the torsion tensor: 



dx^" A dx^T, 



This definition is invariant with respect to general coordinate transformations of x'^ and 
covariant with respect to global rotations. Thus, the torsion tensor has straightforward 
physical interpretation: it is equal to the surface density of the Burgers vector. 
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The physical interpretation of the §0(3) connection is given in section and now we 
show how this definition reduces to the expression for the Burgers vector (jH)) obtained 
within elasticity theory. If the curvature tensor for the SO(3) connection is zero, then, 
according to theorem the connection is locally trivial, and there exists such §0(3) 
rotation such that cj^j-^ = 0. In this case, we return to expression 

If the SO(3) connection is zero and vielbein is a smooth function, then the Burgers 
vector corresponds uniquely to every contour. It can then be expressed as a surface 
integral of the torsion tensor. The surface integral depends only on the boundary contour 
but not on the surface due to the Stokes theorem. 

We have shown that the presence of linear 
defects results in a nontrivial torsion tensor. In 
the geometric theory of defects, the equality of 
the torsion tensor to zero T^jy* = is naturally 
considered the criterion for the absence of dis- 
locations. Then, under the name dislocation 
fall not only linear dislocations but, in fact, ar- 
bitrary defects in elastic media. For example, 
point defects; vacancies and impurities, are also 
dislocations. In the first case we cut out a ball 
from the Euclidean space and then shrink 
the boundary sphere to a point (Fig. Hj). In 
the case of impurity, a point of the Euclidean 
space is blown up to a sphere and the produced 
cavity is filled with the medium. Point defects 
are characterized by the mass of the removed 
or added media, which is also defined by the 
vielbein 
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Figure 4: Point defect: a vacancy ap- 
pears when a ball is cut out from the 
medium, and the boundary sphere is 
shrunk to a point. 
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(10) 



where ?/*(x) are the transition functions to a Cartesian coordinate frame in and po is 
the density of the medium which is supposed to be constant. The mass is defined by the 
difference of two integrals, each of them being separately divergent. The first integral 
equals to the volume of the medium with defects and the second is equal to the volume 
of the Euclidean space. The torsion tensor for a vacancy or impurity is zero everywhere 
except at one point, where it has a 5-function singularity. For point defects, the notion 
of the Burgers vector is absent. 

According to the given definition, the mass of an impurity is positive because the 
matter is added to the media, and the mass of a vacancy is negative because part of 
the medium is removed. The negative sign of the mass causes serious problems for a 
physical interpretation of solutions of the equations of motion or the Schrodinger equation. 
Hence, we make a remark. Strictly speaking, the integral (fTIHl should be called the "bare" 
mass because this expression does not account for elastic stresses arising around a point 
dislocation. The effective mass of such a defect must contain at least two contributions: 
the bare mass and the free energy coming from elastic stresses. The question about the 
sign of the effective mass is not solved yet and demands a separate analysis. 

Surface defects may also exist in three-dimensional space, in addition to point and line 
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dislocations. In the geometric approach, all of them are called dislocations because they 
correspond to a nontrivial torsion. 



4 Disclinations 

In the preceeding section, we related dislocations to a nontrivial torsion tensor. For this, 
we introduced an §0(3) connection. Now we show that the curvature tensor for the SO(3) 
connection defines the surface density of the Frank vector characterizing other well-known 
defects - disclinations in the spin structure of media [HT] . 

Let a unit vector field n*(a;) [n^rii = 1) be given at all points of the medium. For 
example, n* has the meaning of the magnetic moment located at each point of the medium 
for ferromagnets (Fig. Els). For nematic liquid crystals, the unit vector field ra* with the 
equivalence relation ~ — ra* describes the director field (Fig. 



I ^ \ 




Figure 5: Examples of media with the spin structure: (a) ferromagnet, (&) liquid crystal. 

We fix some direction in the medium n\. Then the field rC'ix) at a point x can be 
uniquely defined by the field uj^^x) = — u;-'*(a;) taking values in the rotation algebra so (3) 
(the rotation angle), 

n' = niS/{Lj), 

where 5*/ G SO(3) is the rotation matrix corresponding to the algebra element uj^^ . Here, 
we use the following parameterization of the rotation group §0(3) by elements of its 
algebra (see, e.g., P^) 

S,^ = (e("^))i^' = cosa;5^' + sin a; + ^^(1 - cosa;) G §0(3), (11) 

where {uje)i^ = uj^Ski^ and oj = ^Jufuji is the modulus of the vector u;*. The pseudovector 
Lo^ = where e^^^ is the totally antisymmetric third-rank tensor, e^^'^ = 1, is 

directed along the rotation axis and its length equals the rotation angle. We call the field 
iv^^ spin structure of the media. 

If a media has a spin structure, then it may have defects called disclinations. For 
linear disclinations parallel to the axis, the vector field n lies in the perpendicular plain 
The simplest examples of linear disclinations are shown in Fig. IHl Every linear 
disclination is characterized by the Frank vector 

e, = e,,fc^]^■^ (12) 



10 




Figure 6: The vector field distributions on the plane x^,x'^ for the linear disclinations 
parallel to the axis, (a) |0| = 2tt. (b) |B| = An. 

where 

Q'^ = j dx^d^u'^ (13) 

and the integral is taken along closed contour C surrounding the disclination axis. The 
length of the Frank vector is equal to the total angle of rotation of the field ra* as it goes 
around the disclination. 

The vector field defines a map of the Euclidean space to a sphere n : — >■ S^. For 
linear disclinations parallel to the x^ axis, this map is restricted to a map of the plane 
to a circle In this case, the total rotation angle must obviously be a multiple of 2tt. 

For nematic liquid crystals, we have the equivalence relation ~ — n*. Therefore, for 
linear disclinations parallel to the x^ axis, the director field defines a map of the plane 
into the projective line n : — >■ MP^. In this case, the length of the Frank vector must 
be a multiple of tt. The corresponding examples of disclinations are shown in Fig. [7| 




Figure 7: The director field distribution in the x^,x'^ plane for the linear disclinations 
parallel to the x^ axis, (a) |0| = tt and {b) |0| = Svr. 
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As for the displacement field, the field uj'^^{x), taking values in the algebra so (3), is not 
a smooth function in in the presence of disclinations. We make a cut in M'^ bounded 
by the disclination axis. Then the field uj^^{x) may be considered smooth in the whole 
space except the cut. We assume that all partial derivatives of u!^^{x) have the same limit 
as it approaches the cut from both sides. Then we define the new field 

• • \df,uj'^, outside the cut, 
I lim S^cj*-' , on the cut. 

The functions a;^*-' are smooth everywhere by construction except, may be, on the discli- 
nation axis. Then the Frank vector may be expressed as the surface integral 



= (h dx^'u;^ = dx^ A dx''{d^uj^'^ - d^u;^), (15) 
Jc J J s 

where 5* is an arbitrary surface having the contour C as the boundary. If the field uj^^^ is 
given, then the integrability conditions for the system of equations d^uj'^^ = a;^*-' are 

d^uj^ - d,uj^'^ = 0. (16) 

This noncovariant equality yields the criterion for the absence of disclinations. 

In the geometric theory of defects, we identify the field cj^*-' with the §0(3) connection. 
In the expression for the curvature in ()126|) . the first two terms coincide with (fTBj) . and we 
therefore postulate the covariant criterion of the absence of disclinations as the equality 
of the curvature tensor for the SO(3) connection to zero: 

R>.u'' = 0. 

Simultaneously, we give the physical interpretation of the curvature tensor as the surface 
density of the Frank vector 

(17) 

This definition reduces to the previous expression for the Frank vector ()15|) in the case 
where rotation of the vector n occures in a fixed plane. In this case, rotations are restricted 
by the subgroup §0(2) C SO(3). The quadratic terms in the expression for the curvature 
in (jl26p disappear because the rotation group SO(2) is Abelian, and we obtain the previous 
expression for the Frank vector ()15p. 

Thus, we described the media with dislocations (defects of elastic media) and discli- 
nations (defects in the spin structure) in the framework of Riemann-Cartan geometry, 
the torsion and curvature tensors being identified with the surface density of dislocations 
and disclinations, respectively. The relations between physical and geometrical notions 
are summarized in the Table H 

The same physical interpretation of torsion and curvature was considered in j63j . 
Several possible functionals for the free energy were also considered. In the next section, 
we propose a new expression for the free energy. 
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Existence of defects 






Elastic deformations 








Dislocations 







Disclinations 







Dislocations and disclinations 


7^0 


7^0 



Table 1: The relation between physical and geometrical notions in the geometric theory 
of defects. 



5 Free energy 

Until now, we discussed only the relation between physical and geometrical notions. To 
complete the construction of the geometric theory of defects, we have to postulate equa- 
tions of equilibrium describing static distribution of defects in media. The vielbein e^* 
and §0(3) connection u;^*-' are basic and independent variables in the geometric approach. 
In contrast to previous geometric approaches, we completely abandon the displacement 
field and spin structure as the fields entering the system of equilibrium equations. 
In a general case of a continuous distribution of defects, they simply do not exist. Nev- 
ertheless, at some level and under definite circumstances they can be reconstructed, and 
we discuss this in the following section. 

The expression for the free energy was derived in We assume that equations 

of equilibrium must be covariant under general coordinate transformations and local ro- 
tations, be at most of the second order, and follow from a variational principle. The 
expression for the free energy leading to the equilibrium equations must then be equal to 
a volume integral of the scalar function (the Lagrangian) that is quadratic in torsion and 
curvature tensors. There are three independent invariants quadratic in the torsion tensor 
and three independent invariants quadratic in the curvature tensor in three dimensions. 
It is possible to add the scalar curvature and a "cosmological" constant A. We thus obtain 
a general eight-parameter Lagrangian 

-L = -kR + \Ti,k{PiT''^ + hT^'' + hT^''') 
e 4 

+ \R^MllR''''' + 72^?"''^' + IzR'^'S'') - A, e = det e/, (18) 

where k, /3i,2,3 and 71,2,3 are some constants, and we have introduced the trace of torsion 
tensor Tj = Tij"^ and the Ricci tensor Rik = Rijk' ■ Here and in what follows, transforma- 
tion of the Greek indices into the Latin ones and vice versa is always performed using the 
vielbein and its inverse. For example, 

Rijkl R/iuklf^^ jy TijP; T^ykC^iC j. 

The particular feature of three dimensions is that the full curvature tensor is in a one- 
to-one correspondence with Ricci tensor ()129|) and has three irreducible components. 
Therefore, the Lagrangian contains only three independent invariants quadratic in cur- 
vature tensor. We do not need to add the Hilbert-Einstein Lagrangian R, also yielding 
second-order equations, to the free energy (fTSj) because of identity (jl3(Jj) . 
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Thus, the most general Lagrangian depends on eight constants and leads to very 
complicated equations of equilibrium. At present, we do not know precisely what values 
of the constants describe this or that medium. Therefore, we make physically reasonable 
assumptions to simplify matters. Namely, we require that equations of equilibrium must 
admit the following three types of solutions. 

1. There are solutions describing the media with only dislocations, 
R^J' = 0, T^J ^ 0. 

2. There are solutions describing the media with only disclinations, 
R^J^ + 0, T^,' = 0. 

3. There are solutions describing the media without dislocations and disclinations. 

It turns out that these simple assumptions reduce the number of independent parame- 
ters from eight to two. We now prove this statement. The Lagrangian ()18|) yields the 
equilibrium equations 



1 5L 



e Se^J- 



1 5L 



+ 71 (-^^^,fc^n^^^^■''"e^ + R^''''R^,K^ 
+ 72 [-hl,umR''^'''e\ + R'^'^m,,^^ 
+ 73 (-^i?,fci?^'=e^ + he'R^, + hv^'R^.u^ + Ae^ = 0, 

+ P^-T^e^, + 74v,i?'^^, + l^^.R^r 

+ 73^V, {R\e^, - R\e%) - (z ^ j) = 0, (20) 



(19) 



where the covariant derivative acts with the SO(3) connection on the Latin indices and 
with the Christoffel symbols on the Greek ones. For example. 



■'Vl ^ J"! 

'^Rf'^kj - uj.j'R'"'ik. 



The first condition on the class of solutions of the equilibrium equations is that they 
permit solutions describing the presence of only dislocations in media. This means the 
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existence of solutions with zero curvature tensor corresponding to the absence of disch- 
nations. Substitution of the condition Rfj,y^ = into Eqn for the SO (3) connection 
yields 

(12« + 2/?i-/?2-2/53)T. = 0, 

(/t-/?i-/?2)T* = 0, (21) 
{Ak + 2/5i - (32)Wijk = 0. 

Here T,, T*, and Wijk are the irreducible components of the torsion tensor, 

Tijk = Wijk + T*tijk + -{SikTj — 5jkTi), 

where 

rp^ rp ^ijk rp rp i 

^ '^-'-ijk^ 1 -'- j ij 1 

W,,k = T,,k - T*e,,k - li6^kT, - 6,kT,), W,,ke''' = W,/ = 0. 

In a general case of dislocations all irreducible components of torsion tensor differ from 
zero (Tj, T*, Wijk 7^ 0) and Eqns have a unique solution 

P, = -K, P2 = 2k, /33 = 4«:. (22) 

For these coupling constants, the first four terms in Lagrangian (jl8p are equal to the 
Hilbert-Einstein Lagrangian KR{e) up to a total divergence due to identity (jl3(Jj) . Equa- 
tion (fT^ then reduces to the Einstein equations with a cosmological constant 

RtMu - i^dfivR - l^dfiv = 0. (23) 

In this way, the first condition is satisfied. 

According to the second condition, the equations of equihbrium must allow solutions 
with zero torsion T^^^* = 0. In this case, the curvature tensor has additional symmetry 
Rijki = Rkiij, and Eqn becomes 

(71 + 72 + ^73) V. {R'\e^, - R'^e", - R^"" + i?^^e^) 

+ ^(71 + 72 + 473) {e\e^j - e^e^,) V,R = 0. (24) 

Here, we decompose the Ricci tensor onto irreducible components, 

Rij R ij ~\~ R ~^ ^ '^^^3 ' 



where 



jdS pS p5i ri p^ tdA 

n ij IL ji^ TL i U, Ih ij TL ji. 



Note that for zero torsion, the Ricci tensor is symmetrical: R ij = 0. Contraction of Eqn 
with Cfj,^ leads to the equation 

(71 + 72 + ^73)V,i?^% + ^(71 + 72 + 473)V;,i? = 0. 
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In the general case of nonvanishing curvature, the covariant derivatives V^R^'^ii and V^-R 
differ from zero and are independent. Therefore, we obtain two equations for the couphng 
constants, 

7i + 72 + ^73 = 0, 7i + 72 + 473 = 0. 
which have a unique solution 

7i = -72 = 7, 73 = 0. (25) 

In this case, Eqn (I19|) for the vielbein corresponding to a nonzero torsion also reduces to 
Einstein equations (j2Sl)- 

The last requirement for the existence of solutions with zero curvature and torsion is 
satisfied only for the zero cosmological constant A = 0. 

Therefore, the simple and physically reasonable requirements define the two parameter 
Lagrangian [221 

= -kR + (26) 

which is the sum of the Hilbert-Einstein Lagrangian for the vielbein and the square of the 
antisymmetric part of the Ricci tensor. We note that -R(e) and R^ij{e,u) are constructed 
for different metrical connections. 

Other quadratic Lagrangians were considered, for example, in jnHnBl- We note that 
they also contain the displacement vector as an independent variable along with the 
vielbein. 

Expression defines the free energy density in the geometric theory of defects and 
leads to the equilibrium equations (the Euler-Lagrange equations). In our geometric 
approach, the displacement vector and spin structure do not enter expression for the 
free energy. 

6 Gauge fixing 

In the geometric approach, the vielbein e^* and §0(3) connection cu^j-' are the only vari- 
ables. The displacement field and the spin structure can be introduced only in 
those regions of media where defects are absent. As the consequence of the absence of 
disclinations R^ui-^ = 0, the SO(3) connection is actually a pure gauge ()127|) . i.e., the spin 
structure exists. If, in addition, dislocations are absent (T^u^ = 0) then there is the 
displacement field such that the vielbein equals its partial derivatives ()128j) . In this and 
only in this case can we introduce the displacement field and spin structure. We show 
below that this can be done such that the equations of nonlinear elasticity theory and the 
principal chiral SO(3) field are fulfilled. 

For free energy in ()26p. the Euler-Lagrange equations are covariant under general 
coordinate transformations in M.^ and local SO(3) rotations. This means that any solution 
of the equilibrium equations is defined up to diffeomorphisms and local rotations. For the 
geometric theory of defects to make predictions, we have to fix the coordinate system (to 
fix the gauge in the language of gauge field theory). This allows us to choose a unique 
representative from each class of equivalent solutions. We say afterwards that this solution 
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of the Euler-Lagrange equations describes the distribution of defects in the laboratory 
coordinate system. 

We start with gauge-fixing the diffeomorphisms. For this, we choose the elastic gauge 
proposed in j2S|- This question is of primary importance, and we discuss it in detail. 

Equations (j21), © of elasticity theory for /* = yield the second-order equation for 
the displacement vector 

(1 - 2a)Aui + didju^ = 0, (27) 

where 

A 

2(A + /i) 

is the Poisson ratio, (—1 < cr < 1/2), and A is the Laplace operator. It can be rewritten 
in terms of the induced metric for which we obtain a first-order equation. We choose 
precisely this equation as the gauge condition fixing the diffeomorphisms. We note that 
the gauge condition is not uniquely defined because the induced metric is nonlinear in 
the displacement vector, and different equations for the metric may have the same linear 
approximation. We give two possible choices, 

g'^'^.gup + Y^3^'^^p3,u = 0, (28) 
g^'^y^.gup + Y^^Vp/ = 0, (29) 

where we introduced the notation for the trace of metric = g^^g^y. Gauge conditions 
()28j) and ()29|1 are understood in the following way. The metric g^^, is the Euclidean metric 
written in an arbitrary coordinate system, for example, in the cylindrical or spherical 

o 

coordinate system. The covariant derivative is built from the Christoffel symbols 

corresponding to the metric and V^^^^ = as a consequence. The metric g^y is 
the metric describing dislocations [an exact solution of the equilibrium equations for the 
free energy (j^ ]. The gauge conditions differ because in the first and second cases the 
contraction is performed with the metric of dislocation g'^'^ and the Euclidean metric (7^*^, 
respectively, without changing the linear approximation. Both gauge conditions yield 
Eqn ^TJ\ in the linear approximation in the displacement vector (jSj). This is most easily 
verified in Cartesian coordinates. 

From the geometric standpoint, we have the following. The medium with dislocations 
is diffeomorphic to the Euclidean space equipped with two metrics g^^^ and g^y. The 

metric g^^ is a flat Euclidean metric written in an arbitrary coordinate system. The 
metric g^j,i, is not fiat and describes the distribution of dislocations in the same coordinate 
system. In fact, the metric g^^ is used only to fix the coordinate system in which the 
metric g^j^u is measured. 

If the solution of the equilibrium equations satisfies one of the gauge conditions (j2Hl), 
(j2ni); written, for example, in cylindrical coordinate system then we say that the solution 
is found in the cylindrical coordinates. We suppose here that the distribution of disloca- 
tions in elastic media in the laboratory cylindrical coordinate system is described by this 
particular solution. Analogously, we may seek solutions in a Cartesian, spherical, or any 
other coordinate system. 
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Gauge conditions may also be written for the vielbein e^*, which is defined by Eqn 
()120j) . This involves additional arbitrariness because the vielbein is defined up to local 
rotations. This invariance leads to different linear approximations for the vielbein in terms 
of the displacement vector. We consider two possibilities in Cartesian coordinates: 

e^j ^ - df,Ui, (30) 

ef^i ^ -^{d^Ui + diUf,), (31) 

where the index is lowered with the help of the Kronecker symbol. For these possibilities 
and gauge condition we have two gauge conditions for the vielbein 

g'^'V.e,, + -^e^V^e^ = 0, (32) 
i — Zcr 

g^'V^e^i + -^e'^.V^e^ = 0, (33) 
i — Za 

where = e'^iCf/. These conditions differ in the coefficient before the second term. We 

o 

note that in a curvilinear coordinate system, the covariant derivative must also include 
the flat §0(3) connection acting on indices One can also write other possible gauge 
conditions having the same linear approximation. The question of the correct choice is 
unanswered at present and outside the scope of this review. At the moment, we want 
only to demonstrate that the system of coordinates must be fixed, and that the gauge 
condition depends on the Poisson ratio, which is the experimentally observed quantity. 

Gauge conditions are first order equations by themselves and have some 

arbitrariness. Therefore, to fix a solution uniquely, we must impose additional boundary 
conditions on the vielbein for any given problem. 

If the defect-free case, T^,y* = 0, R/j^uj^ = 0, and the equilibrium equations are satisfied 
because the Euler-Lagrange equations for are satisfied. In this and only in this case, 
we can introduce a displacement vector, and the elastic gauge reduces to the equations of 
nonlinear elasticity theory. In the presence of defects, a displacement field does not exist, 
and the elastic gauge simply defines the vielbein. 

In choosing the free energy functional, we required that the conditions R^u^^ = and 
Tf_iu^ = satisfy the Euler-Lagrange equations. This is important because otherwise we 
would obtain an additional condition on the displacement vector (the Euler-Lagrange 
equations) besides the elasticity theory equations following from the elastic gauge. 

We stress an important point once again. In the geometric theory of defects, we assume 
that there is a preferred laboratory coordinate system in which measurements are made. 
This coordinate system is related to the medium without defects and elastic stresses and 
corresponds to the fiat Euclidean space M^. Gauge conditions (j2HI), and (jS2I), 
are written precisely in this Euclidean space and contain a measurable quantity, the 
Poisson ratio a. This property essentially distinguishes the geometric theory of defects 
from the models of gravity in which all coordinate systems are considered equivalent. 

The elastic gauge is used to fix diffeomorphisms. The expression for the free energy in 
(I26|) is also invariant under local §0(3) rotations, and they must also be fixed. For this, 
we recently proposed the Lorentz gauge for the connection |Mj 

d^tu^/ = 0. (34) 
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This gauge is written in the laboratory Cartesian coordinate system and has deep physical 
meaning. That is, let disclinations be absent {R^uj^ = 0). Then the SO(3) connection is 
a pure gauge: 

u^/ = d^S-'.'^Sk', S/ G §0(3). 
In this case, the Lorentz gauge reduces to the principal chiral SO(3)-field equations 

d^id^S-'/Sk') = 

for the spin structure ci;*-?(x). Principal chiral field models (see, e.g., [64-68]) for different 
groups and in a different number of dimensions attract much interest in mathematical 
physics because they admit solutions of topological soliton types and find broad applica- 
tions in physics. 

Thus, the Lorentz gauge (j34j) means the following. In the absence of disclinations, 
equations of equilibrium are identically satisfied, and there exists a field a;*-' that satisfies 
equations for the principal chiral field. By this we mean that the spin structure of the 
medium is described by the model of the principal chiral field in the defect-free case. 

The principal chiral field model is not the only one that can be used for fixing local 
rotations. The Skyrme model jHSI can also be used for this purpose. The Euler-Lagrange 
equations for this model are not difficult to rewrite in terms of the SO(3) connection and 
use as the gauge conditions. 

There are other models for spin structures. For describing the distribution of magnetic 
moments in ferromagnets or the director field in liquid crystals, one uses the expression 
for the free energy depending on the vector n-field itself jTUl • Lately, much attention 
is paid to the Faddeev model of the n-field jTJ. The question whether there are gauge 
conditions on the SO(3) connection that yield these models in the absence of disclinations 
is unanswered at present. 

Thus, we pose the following problem in the geometric theory of defects: to find the 
solution of the Euler-Lagrange equations for free energy that satisfies the elastic 
gauge for the vielbein and the Lorentz gauge for the §0(3) connection. In sections |H1 
and ^2 we solve this problem for the wedge dislocation in the framework of the classical 
elasticity theory and the geometric theory of defects, respectively, and afterwards compare 
the obtained results. 

7 Asymmetric elasticity theory 

In the preceeding section, we used the elasticity theory and the principal chiral §0(3)-field 
model to fix the invariance of free energy (j^Uj) in the geometric theory of defects. This 
is not a unique possibility, because other models may be used for gauge fixing. In the 
present section, we show how another model - asymmetric elasticity theory - can be used 
for fixing diffeomorphisms and local rotations. 

At the beginning of the last century, the Cosserat brothers developed the theory of 
elastic media, every point of which is characterized not only be its position but also by 
its orientation in space [72], i.e., a vielbein is specified at every point (Fig|H)). 

From physical standpoint, this means that every atom in the crystalline structure is 
not a point but an extended object having orientation. In this case, the stress tensor 
is no longer a symmetric tensor, and the corresponding theory is called the asymmetric 
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elasticity theory. Contemporary exposition of this approach is given in [73J. In the present 
section, we show that the asymmetric theory of elasticity is naturally incorporated into 
the geometric theory of defects. 

The main variables in the asymmetric the- i r 

ory of elasticity are the displacement vector \^ /(^ ▼^X^ *r~~*' 

and the rotation angle lo^'{x). The di- ▼ »■ r - 

rection of the pseudovector cu* coincides with 
the rotation axis of the medium element and 



its length is equal to the angle of rotation. The ▼ Ji ttA ^ 



angle of rotation discussed in section HI is dual 

to the spin structure field u!''^{x): Uij = e^kUJ^. ' 

The Cosserat medium is characterized by ^. o -i^ • <■ i ^ 

, , X 1 r r i'lgure 8: bvery pomt of the Cosserat 

the stress tensor a -'ix) (the density of forces . , . , , , . 

r .,, 1 • • ,1 T medmm is characterized not only by its 

acting on the surface with normal t m the di- ■ ■ ^ i i • ... 

.. .X 1 , , , X position but also by its orientation m 

rection j) and the torque stress tensor ^ ■'[x) 

(the density of torques acting on the surface ^P^'^^- 

with normal i in the direction j). The Cosserat medium is in equilibrium if forces and 
torques are balanced at every point, 

d,a^' + f = 0, (35) 
e'^^ajk + djii^' + = 0, (36) 

where /'(x) and m*(x) are the densities of nonelastic external forces and torques. As 
a consequence of Eqn (j36|l . the stress tensor is symmetric if and only if the condition 
djii^"^ + m* = is satisfied. 

The displacement field and rotation angle uniquely define the deformation tensor ejj(x) 
and the twist tensor Hij{x)\ 

^ij = diUj — Uij, ^^^^ 
l^ij = diUJj. 

In general, the deformation and twist tensors have no symmetry in their indices. 

Hook's law in the Cosserat medium is changed to two linear relations connecting the 
stress and torque tensors with the deformation and twist tensors, 

aij = 2iJ,e{ij} + 2ae[ij] + X^ijek" , (38) 
/iij = 27K{ij} + 2e/t[ij] + pSijKk'', (39) 

where /i, A are the Lame coefficients, and 7,e are four new elastic constants charac- 
terizing the medium. Braces and square brackets denote symmetrization and antisym- 
metrization of indices, respectively. 

In the case where ^ 

^ij = 2^9iUj - djUi), (40) 

the deformation tensor is symmetric and has the previous form (@)). Equation (|36|). to- 
gether with (|HS|) and (jSHI), then reduces to the equation 

(7 + e)e'^''AdjUk + m' = 0. 
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The first term vanishes as a consequence of Eqn (jlTj). Thus, for spin structure 1)401) and 
= we return to the symmetric elasticity theory. 

Equations (jHSj) . ()H^. (jHHj) . and (jSHl), together with the boundary conditions, define 
the equilibrium state of Cosserat media. We now show how this model is included in the 
geometric theory. First, we note that in the absence of defects (T^u^ = 0, Rfj,^^^ = 0), the 
fields and a;*-' exist. Then the vielbein and the SO(3) connection are defined by the 
deformation and twist tensors in the linear approximation: 

e,^ = d.y^S/iu) ^ (5^ - d,u^){6] + uj/) ^ 6; - e,\ (41) 
^ d^uj'^ = e'^'^K.k. (42) 

We note that relations (jHTj) can be regarded as equations for the displacement vector 
and rotation angle. The corresponding integrability conditions were obtained in |74j . 
These integrability conditions are the linear approximations of equalities T^^^* = and 
R^lv^^ = defining the absence of defects. 

If nonelastic forces and torques are absent (/* = 0, m* = 0), then the asymmetric 
theory of elasticity reduces to second-order equations for the displacement vector and 
rotation angle: 

(/i + a) Am' + (/i - a + \)d'djU^ - 2adjU^' = 0, (43) 
(7 + e) Ac^' + (7 - e + p)d'djJ + 2ae'^\d,Uk - uo,k) = 0. (44) 

We rewrite these equations for the vielbein and SO(3) connection 

ifi + a)V%' + ifi-a + A)V'e^ - {fi - a)uj/' = 0, (45) 
^(7 + e)e''''V''u;,,k + \{l-e + f3)e^'^V'u,,k + 2ae'^'e,^ = 0. (46) 

Of course, these are not the only equations that coincide with Eqns (HHjl and (jUj) in 
the linear approximation. At present, we do not have arguments for the unique choice. 
The derived nonlinear equations of the asymmetric elasticity theory can be used as gauge 
conditions in the geometric theory of defects. Here, we have six equations for fixing 
diffeomorphisms (three parameters) and local §0(3) rotations (three parameters). Thus, 
the asymmetric theory of elasticity is naturally embedded in the geometric theory of 
defects. 

In section IHl we considered the elastic gauge for the vielbein and the Lorentz gauge 
for the §0(3) connection. In this case, the spin structure variables do not interact with 
elastic deformations when defects are absent. In the asymmetric elasticity theory, the 
elastic stresses directly influence the spin structure and vice versa. 

8 Wedge dislocation in the elasticity theory 

By wedge dislocation, we understand an elastic medium that is topologically the Eu- 
clidean space without the z = axis - the core of dislocation, obtained in the 
following way. We consider the infinite elastic medium without defects and cut out the 
infinite wedge of angle —2n6. For definiteness, we assume that the edge of the wedge 
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coincides with the z axis (FigE)). The edges of the cut are then moved symmetrically 
one to the other and glued together. After that, the medium moves to the equilibrium 
state under the action of elastic forces. If the wedge is cut out from the medium, then 
the angle is considered negative: — 1 < < 0. For positive ^, the wedge is added. Thus, 
the elastic medium initially occupies a domain greater or lesser than the Euclidean space 
M^, depending on the sign of the deficit angle Q] in the cylindrical coordinates r, yj, z this 
domain is described by the inequahties 



0<r<cxD, 0<(y9< 27ra, — oo < z < oo, 



a 



1 + e. 



(47) 



We note that the wedge dislocation is often 
called disclination. In our approach, this term 
seems unnatural because the wedge dislocation 
is related to a nontrivial torsion. Moreover, the 
term disclination is used for defects in the spin 
structure. 

We now proceed with the mathematical for- 
mulation of the problem for the wedge disloca- 
tion in the framework of elasticity theory. To 
avoid divergent expressions arising for infinite 
media, we suppose that the wedge dislocation 
is represented by the cylinder of a finite radius 
R. This problem has translational symmetry 
along the z axis and rotational symmetry in 
the X, y plane. Therefore, we use cylindrical 
coordinate system. Let 
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Figure 9: Wedge dislocation with the 
deficit angle It^Q. For negative and pos- 
itive 6*, the wedge is cut out or added, 
respectively. 



Ur,U^,U^) " (4Sj 

be the components of the displacement covector with respect to orthonormal basis in the 
cylindrical coordinate system. This covector satisfies the equilibrium equation following 
the substitution of (jHl) in Eqn ^ in domain (jTTjl . 



(49) 



where Vj is the covariant derivative for the flat Euclidean metric in the considered coor- 
dinate system. 

For references, we write expressions for the divergence and Laplacian for a covector 
field in the cylindrical coordinate system: 



AUr 

All,, 



drirdrUr) + —d^^Ur + dlur - -^Ur T^d^U^ 



dr{rdrU^) + 



1 2 2 1 2 



-dr{rdrUz) + -^d Uz + d^u^. 
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Taking the symmetry of the problem into account, we seek the solution of Eqn (jl^ in 
the form 

Ur = u{r), = A{r)ip, Uz = 

where u{r) and A{r) are two unknown functions depending only on the radius. We impose 
the boundary conditions: 

The first four equations are geometrical and correspond to the process of dislocation 
creation. The last condition has simple physical meaning: the absence of external forces 
at the boundary of the medium. The unknown function A{r) is found from the second to 
the last boundary condition (fSUj) 

A{r) = ^ 



1 + 9 



Straightforward substitution then shows that ip and z components of equation of equihb- 
rium are identically satisfied, and the radial component reduces to the equation 

dJrdru) — — = D, D = —- — — — - — = const. 

r 1 — a 1 + 9 

The general solution of this equation is 

D , C2 

u = — r In r + Cir H , Ci 2 = const. 

2 r 

The constant of integration C2 = due to the boundary condition at zero. The constant Ci 
is found from the last boundary condition in Finally, we obtain the known solution 
of the considered problem [^2] 

D , r 

2 ^ (51) 

The letter e in the expression for Ur denotes the base of the natural logarithm. We note 
that the radial component of the displacement vector diverges as R ^ 00. This means 
that the description of the wedge dislocation requires considering a finite-radius cylinder. 

The linear elasticity theory is applicable for small relative displacements, which for a 
wedge dislocation are equal to 

dur 9 1 — 2cr , r 1 du^ 9 

— In — 



dr 1 + ^2(1 - a) R r dip 1 + 9 

This means that we are able to expect correct results for the displacement field for small 
deficit angles {9 <^ 1) and near the boundary of the cylinder (r ~ i?). 

We find the metric induced by the wedge dislocation in the linear approximation in 
the deficit angle 9. Calculations can be performed using the general formulas (0) or the 
known expression for the variation of the form of the metric 

o o 

Sgf^iy = -Vp^U^ - VyU^. (52) 
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After simple calculations, we obtain the expression for the metric in the x, y plane: 

ds'=(l + e^-f^ In ^ dT' + T'(l + e^-f^ In ^ + e-^\ V- (53) 
y 1 — (T Rj \ 1 — a R \ — o) 

This metric is compared with the metric obtained as the solution of three-dimensional 
Einstein equations in section ITTl 

9 Edge dislocation in the elasticity theory 

Wedge dislocations are relatively rarely met in nature because they require a much amount 
of a medium to be added or removed, resulting in a vast quantity of energy expenses. 
Nevertheless, their study is of great importance because other linear dislocations can be 
expressed as a superposition of wedge dislocations. In this respect, wedge dislocations 
are elementary. We show this for an edge dislocation - one of the most widely spread 
dislocations, as an example. An edge dislocation with the core coinciding with the z axis 
is shown in FigJTUla. It appears as the result of cutting the medium over the half-plane 




Figure 10: Edge dislocation with the Burgers vector h directed to the dislocation axis 
(a). Edge dislocation as the dipole of two wedge dislocations with positive and negative 
deficit angles (6). 

2/ = 0, X > 0, moving the lower edge of the cut towards the z axis on a constant (far from 
the core of the dislocation) Burgers vector b, and subsequently gluing the edges. To find 
the displacement vector field for the edge dislocation, we may solve the corresponding 
boundary value problem for the equilibrium equations pH)) [61 . However, we follow 
another way, knowing the explicit form of the displacement vector for a wedge dislocation. 
The edge dislocation is represented by the dipole of two wedge dislocations with positive, 
27r^, and negative, —2ti6, deficit angles as shown in FigllOb. We assume that the axes of 
the first and second wedge dislocations are parallel to the z axis and intersect the x, y plane 
at points with the respective coordinates (0, a) and (0, — a). The distance between the 
wedge dislocation axes is equal to 2a. It follows from the expression for the displacement 
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field 1)511) that far away from the origin (r ^ a), the displacement field for the wedge 
dislocations has the following form in the first order in small 9 and a/r: 
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It is sufficient to sum displacement fields (fS^ and to find the displacement field 
for the edge dislocation because the elasticity theory equations are linear. After simple 
calculations, up to the translation of the whole medium by a constant vector along the y 



axis, we obtain 



Uy = —b 
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(56) 



2(1 -a) eR 2{1 - a) x^ + _ 
where we have introduced the notation for the modulus of the Burgers vector 

b=\b\ = -2a9. 

This result coincides with the expression for the displacement field obtained by direct solu- 
tion of the elasticity theory equations [HI]. Thus, we have shown that an edge dislocation 
is the dipole of two parallel wedge dislocations. 

We next find the metric induced by the edge dislocation. Using formulas ()52|). we 
obtain the metric in the x,y plane in the linear approximation in 6 and a/r: 



1 — 2a 6 . \ / , 2 1 1 i\ 26 cos 09 , , 
1 + sin [dr + r dip ) ^ —drd^p. 



1 — a r J 1 — a 

We note that the induced metric for an edge dislocation does not depend on R. 



(57) 



10 Parallel wedge dislocations 

In the absence of disclinations {Rf_iui'' = 0), the SO(3) connection is a pure gauge, and 
equations of equilibrium for the SO(3) connection are identically satisfied. In this 
case, the explicit form of the §0(3) connection is uniquely determined by the spin struc- 
ture u"^^ . The field uj^^ satisfies equations for the principal chiral field as a consequence of 
the Lorentz gauge in (jM)). Solution of this system of equations defines the trivial SO(3) 
connection. Thus in the absence of disclinations, the problem is reduced to solution of 
the Einstein equations for the vielbein in the elastic gauge and solution of the principal 
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chiral field model for the spin structure. After that, we can compute the torsion tensor 
through the formulas ()125|) . which defines the surface density of the Burgers vector. 

Because disclinations are absent (the curvature tensor is equal to zero), we have the 
space of absolute parallelism. The whole geometry is then defined by the vielbein e^', 
which uniquely defines the torsion tensor via (jl25|) for a vanishing §0(3) connection. Here, 
we assume that a trivial SO(3) connection is equal to zero. The vielbein e^* satisfies the 
three-dimensional Einstein equations with a Euclidean signature metric, which follow from 
the expression for the free energy for R^uj'' = 0, 

~ 1 

Rfiu — -^9iJ.vR = Tfj,y. (58) 

Here, we have added the source of dislocations T^j^ to the right hand side of Einstein 
equations (it is the energy- momentum tensor in gravity). 

We note that without a source of dislocations, the model would be trivial. Indeed, 
the scalar curvature and Ricci tensor are equal to zero: i? = 0, R^j^u = for T^i, = 
as a consequence of Einstein equations (|HH|l . Then, the full curvature tensor without 
sources is also equal to zero because in three-dimensional space, it is in a one-to-one 
correspondence with Ricci tensor p29|) . The vanishing of the full curvature tensor means 
the triviality of the model because defects are absent in this case. The similar statement in 
three-dimensional gravity is well known. It is usually formulated as: "Three-dimensional 
gravity does not describe dynamical, i.e., propagating degrees of freedom". 

For our purposes, we have to find the solution of Einstein equations ()58|) describing 
one wedge dislocation. The Einstein equations are a system of nonlinear second-order 
partial differential equations. Not too many exact solutions are known at present, even in 
a three-dimensional space. The remarkable exact solution describing an arbitrary static 
distribution of point particles is well known in three-dimensional gravity for the Lorentz 

signature metric (H ) [75-77]. We find this solution for the Euclidean signature metric 

and show that it describes an arbitrary distribution of parallel wedge dislocations in the 
geometric theory of defects. Hence, we first consider the more general case of an arbitrary 
number of wedge dislocations and then analyze in detail one wedge dislocation which is 
of interest to us. We do this deliberately because the solution in a more general case 
does not involve essential complications. At the same time, an arbitrary distribution 
of wedge dislocations is much more interesting for applications. For example, the edge 
dislocation was shown in the preceeding section to be represented by a dipole of two wedge 
dislocations of different signs. 

We consider elastic medium with arbitrary distributed but parallel wedge dislocations. 
We choose the coordinate system such that the z = axis is parallel to dislocations axes, 
and axes = {x,y}, a = 1,2 are perpendicular to the z axis. Then the metric has 
the block-diagonal form 

ds"^ = + N'^dz^, (59) 

where 

dt^ = gafjdx^'dx^ 

is a two-dimensional metric on the x,y plane. A two-dimensional metric gaj3{x,y) and a 
function N{x, y) are independent of z due to translational symmetry along the z axis. 

We can say this more simpler dropping the physical arguments. We consider the 
block-diagonal metric of form (j59j) . which has translational invariance along the z axis. 
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Afterwards, we show that the corresponding solution of the Einstein equations indeed 
describes an arbitrary distribution of parallel wedge dislocations. 
The curvature tensor for metric ()59|) has the components: 



p Z p 5 ri 



1 



(2) r 

where R^p^ is the curvature tensor for the two-dimensional metric gap and Vq is the 
two-dimensional covariant derivative with the Christoffel symbols. 
We choose the source of dislocations as 

9-rr ^ 

(60) 

Tafi Tpfz 0, 

where 5*^^)(r — r„) = 5{x — Xn)5{y — Un) is the two-dimensional (5-function on the x^y plane 
with the support at a point r„ = {xn,yn)- The factor g^"^^ = del gap in front of the sum 
sign is due to the property of the 5-function, which is not a function but a tensor density 
with respect to general coordinate transformations. We show later that the solution of 
the Einstein equations with such a source describes M parallel wedge dislocations with 
deficit angles 9n, which intersect the x, y plane at the points yn)- In three-dimensional 
gravity this source corresponds to particles of masses nin = 27r9n, being at rest at the 
points r„. 

Einstein equations (fSHj) then reduce to four equations, 

V„V^iV - gafsV^V^N = 0, (61) 

-iiV^/?(^) = ^5^M(^)(r-r„), (62) 

V 9 n=l 

where R^"^^ is the two-dimensional scalar curvature. 

The metric of form (j59p is still invariant under coordinate transformations in the x, y 
plane. Using this residual symmetry, we fix the conformal gauge on the plane (this is 
always possible locally) 

gafS = e^'^dap, 

where (f){x,y) is some function. 

In the conformal gauge, Eqn ()61|) becomes 

dadfsN = 0. 

For constant boundary conditions for on the boundary of the x, y plane this equation 
has the unique solution = const. Changing the scale of z coordinate, we can set A^ = 1 
without loss of generality. Then Eqn (jH^ reduces to the Poisson equation 
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which has the general solution 



^ In Ir — r„| + ^ In C, C = const > 0. 



2 

n 

Thus the metric in the x, y plane is 

(iZ^ = CjJ|r-r„p^"(dr2 + r2d(^2^, < r < oo, < < 27r, (63) 

n 

where the polar coordinates r, ip cover the whole plane and not more (this is important 
!). Any solution of the Einstein equations is defined up to choosing the coordinate system 
because the equations are covariant. Using this, we set C = 1, which is always possible 
by choosing the scale of r. 

This is indeed the exact solution of the Einstein equations describing arbitrary dis- 
tribution of parallel wedge dislocations. This statement is made clear from the following 
consideration. 

We note that transition to a continuous distribution of dislocations in the geometric 
approach is simple. For this, we have to substitute a continuous distribution of sources 
in the right hand side of the Einstein equations instead of 5-sources. 

We consider one wedge dislocation with the source at the origin in more detail to show 
that metric ()63|) indeed describes an arbitrary distribution of wedge dislocations 

T^. = ^=e5^'\x.y). (64) 

The corresponding metric (j63|) for C = 1 is 

dl^ = r^\dr^ + r^dip^). (65) 
We pass to a new coordinate system 

r' = -r", if' = aifi, a = 1 + 9, (66) 
a 

in which the metric becomes Euclidean 

df = dr'^ + r'^d^'\ (67) 

but the range of the polar angle differs now from 27r: < < 27rQ;, and covers the x,y 
plane with removed or added angle 2tt9. 

Because the metric coincides with the Euclidean one in the primed coordinate system 
r', ip' , we have the Euclidean plane with a removed or added wedge because the angle ip' 
varies within the interval (0,27ra). The transformation to coordinates r, in ()66|1 means 
the gluing of the edges of the wedge that has appeared, which produces a cone. Therefore, 
both metrics and ()67p describe the same geometric object - conical singularity. The 
torsion and curvature tensors are obviously zero everywhere except at the origin. 

Creation of a conical singularity coincides exactly with creation of the wedge disloca- 
tion in the geometric theory of defects. It is not difficult to show that general solution (j63j) 
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describes an arbitrary distribution of conical singularities with deficit angles 6.^ at points 
Vn. Thus, this solution describes an arbitrary distribution of parallel wedge dislocations. 

In the next section, we consider one wedge dislocation in detail. For this, we perform 
one more coordinate transformation 

/ = ar', ^ = -if'. (68) 

a 

Metric fl67|) then becomes 

df = ^df + fd^^, a = 1 + 6. (69) 
This is one more frequently used form of the metric for a conical singularity. 



11 Wedge dislocation in the geometric approach 

We now consider a wedge dislocation from the geometric standpoint. From the qualitative 
standpoint, the creation of a wedge dislocation coincides with the definition of conical 
singularity. However, there is a quantitative difference because metric fl69j) depends only 
on the deficit angle 6 and cannot coincide with the induced metric obtained within 
the elasticity theory. The difference arises because we require the displacement vector 
in the equilibrium to satisfy equilibrium equations after removing the wedge and gluing 
its edges (creating conical singularity) in the elasticity theory. At the same time the 
X, y plane for a conical singularity after the gluing can be deformed in an arbitrary way. 
Formally, this manifests itself in that metric (j53|) obtained within the elasticity depends 
explicitly on the Poisson ratio, which is absent in gravity theory. 

We proposed the elastic gauge for solving this problem ^23 . We choose elastic gauge 
()33p as the simplest one for a wedge dislocation. This problem can be solved in two ways. 
First, the gauge condition can be inserted in the Einstein equations directly. Second, 
we can find the solution in any suitable coordinate system and then find the coordinate 
transformation such that the gauge condition is satisfied. 

It is easier to follow the second way because the solution for the metric is known, Eqn 
(jUUj) . The vielbein can be associated with metric (jU^ as 

a 

Here, a hat over an index means that it corresponds to the orthonormal coordinate system, 
and an index without a hat is the coordinate one. Components of this vielbein are the 
square roots of the corresponding metric components and therefore have the symmetric 
linear approximation ()31|) . We transform the radial coordinate / — > /(r) because a wedge 
dislocation is symmetric under rotations in the x, y plane. After that transformation the 
vielbein components take the form 

e/ = -, e/ = /, (70) 
a 

where the prime denotes differentiation with respect to r. We choose the vielbein corre- 
sponding to the Euchdean metric as 

e/ = 1, e/ = r. (71) 
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It defines the Christoffel symbols F^j^'' and SO(3) connection u^i^ , which define the co- 
variant derivative. We write only nontrivial components 

o o 1 o 

-p v> — -p V> — _ r r _ 

° _ _ ° J — ^ 

Substitution of the vielbein into gauge condition yields the Euler differential equation 
for the transition function 

a ar 1 — 2cr \ a r r^y 

Its general solution depends on two arbitrary constants Ci^2, 

/ = C'ir^i+C2r^2, (72) 
where the exponents 71^2 are defined by the quadratic equation 

7^ + 27^6 - a = 0, b ^ 



2(1 -a) 



which has real roots for 6 > —1 with different signs: positive root 71 and negative 72. 
We recall that there are thermodynamical constraints — l<cr<l/2on the Poisson ratio 

To fix the constants, we impose boundary conditions on the vielbein: 

^r\=n = h e/l^^ = 0. (73) 

The first boundary condition corresponds to the last boundary condition for displacement 
vector (jSUj) (the absence of external forces on the surface of the cylinder), and the second 
one corresponds to the absence of the angular component of the deformation tensor at 
the core of dislocation. Equations (f7H|) define the values of integration constants 

Ci = — C2 = 0. (74) 

The obtained vielbein defines the metric 

df = ( ^Y'"^' ( dr' + , (75) 



where 



7i 



-eb + veW+T+e. 



This is the solution of the posed problem. The derived solution is valid for all deficit 
angles 6 and for all < r < R. The obtained metric depends on three constants: 6, a, 
and R. The dependence on the deficit angle 6 is due to its occurrence on the right hand 
side of Einstein equations (j^Hj) . The dependence on the Poisson ratio comes from elastic 
gauge (jnni), and, finally, the dependence on the cylinder radius comes from boundary 
condition (|7H|l . 
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If a wedge dislocation is absent, then 6' = 0, a = 1, 71 = 1, and metric (f73j) goes to 
the Euchdean one dl"^ = dr"^ + r^dip"^, as expected. 

We compare metric (f75j) obtained within the geometric approach with the induced 
metric from elasticity theory in Eqn (j53|) . First, it has a simpler form. Second, in the 
linear approximation in 6 

and metric (j75p can be easily shown to coincide precisely with metric (|53|) obtained 
within the elasticity theory. We see that induced metric (j53|) provides only the linear 
approximation for the metric obtained within the geometric theory of defects, which, in 
addition, has a simpler form. Beyond the perturbation theory, we see essential differences. 
In particular, metric (j53|l is singular at the origin, whereas metric (|T5|l obtained beyound 
the perturbation theory is regular. 

The stress and deformation tensors are related by Hook's law 0. There is an ex- 
perimental possibility to check formulas ()75|) because the deformation tensor is the linear 
approximation for the induced metric. For this one has to measure the stress field for a sin- 
gle wedge dislocation. In this way, the geometric theory of defects can be experimentally 
confirmed or discarded. 

The problem of reconstruction of the displacement field for a given metric reduces to 
solution of differential equations (jSj) with metric ()75|) on the right hand side and boundary 
conditions (jSOl)- We do not discuss this problem here. We note that in the geometric 
theory of defects, a complicated stage of finding the displacement vector where it exists 
is simply absent and is not necessary. 

Two-dimensional metric (f7H|l describes the conical singularity in the elastic gauge. The 
relation between conical singularities and wedge dislocations was established in [22,78-80]. 
In these papers, the metric was used in the other gauges (coordinate systems). 



12 Elastic oscillations in media with dislocations 

Elastic oscillations in elastic media without defects are described by a time-dependent 
vector field u\t,x) that satisfies the wave equation (see, e.g., |ST] ) 

Pou' - pAu' - (A + p)d'djU^ = 0, (76) 

where dots denote differentiation with respect to time and po is the mass density of 
medium, which is assumed to be constant. If the medium contains defects, then the 
metric of the space becomes nontrivial, 6ij g^^ = e^^ejdij. We assume that relative 
displacements for elastic oscillations are much smaller than stresses induced by defects: 

d,u' « e/. (77) 

Then, in the first approximation, we assume that elastic oscillations propagate in a Rie- 
mannian space with a nontrivial metric induced by dislocations. Here, we discard changes 
in the metric due to elastic oscillation themselves. Therefore, for elastic oscillations we 
postulate the following equation, which is a covariant generalization of (f7B|) for spatial 
variables: 

Pqu' - fiAu' - (A + p)V'VjU^ = 0, (78) 
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where denote components of displacement vector field with respect to an orthonormal 
space basis Cj (see the Appendix), A = V*Vj is the covariant Laplace-Beltrami operator 
built for the vielbein e^*, and Vj is the covariant derivative. The explicit form of the 
covariant derivative of the displacement field is 

where uj^^^ is the SO(3) connection built for zero torsion. 

We note that the displacement vector field describing elastic oscillations is not the 
total displacement vector field of points of a medium with dislocations. It was already 
said in section El that the displacement field for dislocations can be introduced only in 
those regions of media where defects are absent. If we denote it by mJ,, then the total 
displacement field in these regions is defined by the sum + -u*. There, the vielbein 
is defined only by the displacement field for dislocations e^* = d^u\^. We note that the 
smallness of relative deformations (f77j) is also meaningful in those regions of space where 
displacements -u^ are not defined. 

We now decompose the displacement field covariantly into transversal u^^ and longi- 
tudinal parts m", 

e = u"' + m", 

which are defined by the relations 

V,n"^ = 0, (79) 
Viu) - Vju\ = 0. (80) 

The decomposition of a vector field into longitudinal and transversal parts in three- 
dimensional space with an accuracy up a constant is unique. We recall that Latin indices 
are lowered with the help of Kronecker symbols, = Wj, and this operation commutes 
with covariant differentiation. Equation (j8(Jj) can be rewritten as 

Viu) - Vju\ = e>'ie''j{V^ul - V.u);) = e^ic" j{d^ul - d.u';,) = 0, 

because transformation of Latin indices into Greek ones commutes with the covariant 
differentiation, and the Christoffel symbols are symmetrical in the first two indices. The 
last equality means that the 1-form dx'^u^^ is closed. It is easily confirmed that Eqn 
(f7H|) for elastic oscillations is equivalent to two independent equations for transverse and 
longitudinal oscillations, 

where 

2 ^ _fJ_ 2 ^ \ + 2fl 

Po Po 

are the squares of sound velocities for transverse and longitudinal oscillations. 

Particles arising after the secondary quantization of Eqns (jHTj) are called phonons in 
solids. Therefore, strictly speaking, the problem of scattering of phonons on dislocations 
is a quantum mechanical one. In the present review, we consider only classical aspects of 
this problem. 
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Wave equations (jHTj) contain second and first derivatives of tlie displacement field. 
The latter are contained in the covariant Laplace-Beltrami operator A. The terms with 
second derivatives can be written in the four- dimensional form g'^^dadj^ where g'^^ is the 
inverse metric to 

^ (o -L) • («^) 

where c is either the transverse or longitudinal sound velocity. Above, we used the follow- 
ing notations. Four-dimensional coordinates are denoted by Greek letters from the be- 
ginning of the alphabet {a;°} = = t, and letters from the middle of Greek 
alphabet denote only space coordinates {x'^} = {x^,x'^,x^}. This rule can be easily re- 
membered by the following inclusions {1, 2, 3} C {0, 1, 2, 3} and {/i, u, . . .} C {a, /?,...}. 
Christoffel symbols can be computed for four-dimensional metric ()82j) . which defines 
a system of ordinary nonlinear equations for extremals x"(r) (lines of extremal lengths 
that coincide with geodesies in Riemannian geometry), where dots denote differentia- 
tion with respect to canonical parameter r. For the block-diagonal metric in ()82|) . these 
equations decompose: 

x° = 0, (83) 
x" = -f^p^x'^xP, (84) 

where T^p'^ are the three-dimensional Christoffel symbols constructed for the three-dimensional 
metric g^j^i, which, as we recall, depends only on spatial coordinates for a static distribu- 
tion of defects. Let {x°'{t)} be an arbitrary extremal in four-dimensional space-time. For 
metric (jH2I), its natural projection on space {x"(r)} {0,x^^{t)} is also an extremal but 
now for the spatial part of the metric g^^iy. 

Equations for extremals (jHHj) . (jH^ are invariant under linear transformations of the 
canonical parameter r. Therefore, the canonical parameter can be identified with time, 
T = t = x^, without loss of generality as a consequence of Eqn (jH^ . 

We assume that a particle moves in space along an extremal x^{t) with velocity v. 
This means that 

g ^f^^ ^ ^ • 

The length of the tangent vector to the corresponding extremal {t, x^{t)} in four-dimensional 
space-time is then equal to 

Hence, if the particle moves in space along extremal with velocity less than, equal to, 
or greater than the speed of sound, then its world line in space-time is timelike, null, or 
spacelike, respectively. 

We return now to propagation of phonons in media with defects. As in geometric optics 
[HI], there are useful notions of wave fronts and rays in the analysis of the asymptotic 
form of solutions for wave equations ()81|) . We do not consider mathematical aspects of 
this approach, which is nontrivial and complicated JH^, and instead give only a physical 
description. In the eikonal (high frequency) approximation, phonons propagate along 
rays coinciding with null extremals for the four- dimensional metric ga^^. Forms of rays, 
which are identified with trajectories of phonons, are defined by the three-dimensional 
metric g^^,. This means that in the eikonal approximation, trajectories of transverse and 



33 



longitudinal phonons in a medium with defects are the same and are defined by Eqn ()84|) . 
The difference reduces to the velocities of propagation for transverse and longitudinal 
phonons being different and equal to Ct and Cl, respectively. 



13 Scattering of phonons on a wedge dislocation 

Calculations in the present section coincide, in fact, with the analysis performed in Section 
3 of jSni- The difference is that in what follows, we use the metric written in the elastic 
gauge. This is important because we assume that trajectories of phonons seen in an 
experiment coincide with extremals for the metric precisely in this gauge. 

In the presence of one wedge dislocation, the space metric in the cylindrical coordinates 
r, 99, z is 
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where the nontrivial part of the metric in the r, plane was obtain earlier [see Eqn (j75|) ]. 
Here, we change 71 to 7 for simplicity of notation. The inverse metric is also diagonal. 
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Christoffel symbols for metric ()85j) are calculated according to formulas As a 

result, only four Christoffel symbols differ from zero: 



r 
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7-1 



2 

a r 



V = r 
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r 



In the preceeding section, we showed that in the eikonal approximation, phonons 
propagate along extremals x^(t) defined by Eqns (fH^ . In the case considered here, these 
equations are 



7 - 1,2 , -2 



-r H np 

r 7 



27 . . 

if = rip, 

r 

z = 0, 



(86) 
(87) 



where the dot denotes differentiation with respect to time t. It follows from the last 
equation that phonons move along the z axis with constant velocity, which corresponds 
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to translational invariance along z. This means that scattering on a wedge dislocation is 
reduced to a two-dimensional problem in the r, plane, as could be expected. 

The system of equations for r(t) and </9(t) in ()86|) and ()87|) can be explicitly integrated. 
For this, we find two first integrals. First, for any metric, there is the integral for equations 
for extremals 

g^ivX^x^ = const. 

We then have the equality 

2 

^27-2^2 ^ " ^27^2 = Cq = COHSt > 0. (89) 

7 

Second, the invariance of the metric under rotations about the z axis results in the 
existence of an additional integral. It is constructed as follows. There is a Killing vector 
corresponding to the invariance of the metric, which in cylindrical coordinates has the 
simple form k = d^p. Straightforward verification proves that 

g^uk^i''^ = const. 

In the considered case, this results in the identity 

r2T0 = Ci= const. (90) 

We analyze the form of an extremal r = r (</?). First derivatives can be found from 
Eqns dHHl) and P^l 



±r-2^+^/Cor27-^C?, (91) 
= Cir-^\ (92) 

Admissible values of the radial coordinate r for which the expression under the square root 
is nonnegative are to be found later. From the above equations, we obtain the equation 
defining the form of nonradial (Ci 7^ 0) extremals 



This equation can be easily integrated, and we finally obtain explicit formulas defining 
the form of an extremal: 



— s\n^[a{ip + ifo)] = I, (94) 
r 



where 

1/7 

= const > 0, v^o = const. 



( C,a Y 
Vv^7/ 



The constant is positive and defines the minimal distance at which an extremal ap- 
proaches the core of the dislocation, i.e., r > r^. Only for these values of r is the 
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expression under the root in Eqn ()93|) nonnegative. The integration constant (fo is arbi- 
trary and corresponds to the invariance of the problem under rotations around the core 
of dislocation. 

Equations for extremals (jHU)) and (|H7|) also have degenerate solutions: 

—r'^ = ±^/C^{t + to), (fi = const, to = const. (95) 

7 

These extremals correspond to the radial motion of phonons. Such trajectories are un- 
stable in the sense that there are no nonradial extremals near them. 

We note that circular extremals r = const are absent as a consequence of Eqn (j86|) . 
although integrals of motion (j89|) and (j9(J|) admit such a solution. This is because Eqn 
()86p was multiplied by r in deriving first integral (jHH)- 

We now analyze the form of nonradial extremals ()94|) . For any extremal, the radius r 
decreases first from infinity to the minimal value and then increases from to infinity. 
We can assume here without loss of generality that the argument of sine in (j94|) varies 
from to vr. Thus, we obtain the range of changes of the polar angle 

TT 

< ip + ipo < —. 

a 

This means that the extremal comes from infinity at the angle —ipo and goes to infinity 
at the angle n/a — ipo. This corresponds to the scattering angle 

X = vr - - = — 96 
a 1 + o 

We note that the scattering angle depends only on the deficit angle 9 and does not de- 
pend on the elastic properties of the medium. The scattering angle has a simple physical 
interpretation. For positive 6, the medium is cut and moved apart. A wedge of the same 
medium without elastic stresses is inserted in the cavity that has appeared. Afterwards, 
gluing is performed, and the wedge is compressed. The compression coefficient for all cir- 
cles centered at the origin is equal to 1/(1 + 6') due to symmetry considerations. Therefore, 
the scattering angle equals half the deficit angle times the compression coefficient, 

2n9 1 
X = ^ X 



1 + 9 

For ^ = 0, the dislocation is absent, and the scattering angle is equal to zero. 

For positive deficit angles, the scattering angle is positive, which corresponds to re- 
pulsion of phonons from the wedge dislocation. Corresponding extremals are shown in 
FiglTTb. and they have asymptotes as r — > oo. We note that for positive deficit angles, no 
two points on the r, (p plane can be connected by an extremal, i.e., there is a domain to 
the right of the wedge dislocation that cannot be reached at all by a phonon falling from 
the left. 

All extremals shown in the figures in the present section are calculated numerically. 
The values of the deficit angle 9, the scattering angle Xy ^ind the minimal distance 
to the dislocation axis are shown in the figures. For definiteness, we choose cpQ = n 
corresponding to the fall of phonons from the left. 
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Figure 11: (a) Extremals for the wedge dislocation with a positive {6 > 0) deficit angle. 
Two extremals and their reflections with respect to the x axis are shown for the same ^ > 
but different r^- (b) Extremals for the wedge dislocation with a negative (—1/2 < 6 < 0) 
deficit angle. Two extremals are shown for the same but different deficit angles 9. 

For negative deficit angles the scattering angle is negative and is defined by the same 
formulas (jUBj) . This corresponds to the attraction of phonons to the dislocation axis. In 
Figlllb. we show two extremals with the same parameter but for two dislocations 
with different deficit angles. For —1/2 < 6 < 0, the scattering angle varies from to 27r 
(Figlllb). For 6 = —1/2, the scattering angle is equal to 27r. We note that for negative 
deficit angles, phonons have no asymptotes as r — > oo, i.e., for 6 = —1/2, phonons fall 
from infinity {x —oo,y +oo) and return to the infinity {x —oo,y —oo). 

When the deficit angle is sufficiently small (—1 < 6 < —1/2), a phonon makes one 
or several rotations around the dislocation and then goes to infinity. Examples of such 
trajectories are shown in Figs. 11^141 
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Figure 12: For 6 = —2/3, an extremal makes one rotation around the dislocation and 
then goes forward in the original direction. The blown-up part of the trajectory is shown 
in the square to the right. 




Figure 13: For 9 = —3/4, an extremal makes two turns around the dislocation and then 
goes back. The blown-up part of the trajectory is shown in the square to the right. 
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Figure 14: For 9 = —4/5, an extremal makes two and a half rotations around the 
dislocation and then goes forward in the original direction. Subsequent blown-up parts 
of the trajectory are shown in the two squares to the right. 

We consider the asymptotic behavior of nonradial extremals as r oo. As a conse- 
quence of Eqn (PT|) . far from the core of dislocation, we obtain 



-7+1 



It follows from this equation that the dependence of the radius on time is the same as for 
radial extremals . Because 7 > 0, an infinite value of r corresponds to an infinite value 
of time t. This means that the r, (p plane with the given metric is complete at r — 00. 
The origin (the core of dislocation) is a singular point. Only radial extremals fall into it 
at a finite moment of time. 

Integrals of motion (j89|l and (j9n|) have simple physical meaning. Equations for ex- 
tremals (|84p follow from the variational principle for the Lagrangian 

L = ^g^.x'^x", (97) 

describing motion of a free massless point particle in a nontrivial metric gp^p{x). Here, the 
metric is considered as a given external field and is not varied. 
The energy corresponding to this integral is equal to 

1 -u-u 1 -2 -2 1 -2 

If the metric, as in our case, does not depend on time explicitly, then the energy is con- 
served {E = const) and its numerical value for the motion in the r, ip plane is proportional 
to the integral of motion Co- 

For a wedge dislocation, the metric is independent on the polar angle if, and the 
Lagrangian is invariant under rotations: ip ip + const. By the Noether theorem, the 
angular momentum conservation law corresponds to this invariance. 
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As a consequence, the constant of integration Ci is proportional to the angular momen- 
tum. 

We note that the behavior of extremals differs qualitatively from trajectories of point 
particles moving in flat space with Euclidean metric 6^^, in an external potential field U (x) . 
It can be easily shown that trajectories of point particles described by the Lagrangian 

L = —Sfj^yX^x U 

cannot coincide with extremals (j5^ for any function U{x). 

In the present section, we have showed how the problem of scattering of phonons 
on the simplest wedge dislocation is solved in the geometric approach. The results of 
calculations can be checked experimentally. The more general problem of scattering of 
phonons on an arbitrary distribution of wedge dislocations, including the edge dislocation 
and continuous distribution of dislocations, is solved in jSH]. Solution of this problem in the 
geometric approach is simpler than in the framework of classical elasticity theory, where 
we have to solve partial differential equations with complicated boundary conditions. In 
[HH] . calculations were performed in the conformal gauge for the metric. For comparison 
with experiment, the results must be recalculated in the elastic gauge, which has physical 
meaning in the geometric approach. 



14 An impurity in the field of a wedge dislocation 

We consider elastic media with one wedge dislocation containing one atom of impurity 
or vacancy. If we consider the influence of impurity on the distribution of elastic stresses 
small compared with the elastic stresses induced by the dislocation itself, then the motion 
of the impurity can be considered as taking place in three-dimensional space with the 
nontrivial metric (f73jl . In the geometric approach, we assume that the potential energy of 
the interaction equals zero, and all interactions are due to changes in the kinetic energy, 
which depends explicitly on the nontrivial metric. 

We solve the corresponding quantum mechanical problem. We consider bounded states 
of impurity moving inside a cylinder of radius R in the presence of a wedge dislocation. 
We assume that the cylinder axis coincides with the core of dislocation. The stationary 
Schrodinger equation is 

= E^, (98) 

where h is the Plank constant and M, \E', and E are mass, wave function, and energy of the 
impurity. The nontriviality of the interaction with dislocation reduces to the nontrivial 
Laplace-Beltrami operator 

where the metric was found earlier [see Eqn (ffH|) ]. and g = detg^^. 

Taking the symmetry of the problem into account, we solve the Schrodinger equation 
in cylindrical coordinates by the separation of variables. Let 



^{r,ip,z) = Z{z) J2 ^m{r)e 



inup 
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where we have the two following possibilities for the normalization function Z{z). If the 
impurity moves freely along the z axis with momentum hk, then 

Z(z) = e'^'. 



If its motion is restricted by the planes ^ = and z = Zq, then 

Z{z) = J— sin {kiz), ki = —. 

V Zo Zq 

In what follows, we drop the integer-valued subscript / indicating the restricted motion, 
having both possibilities in mind. 

The condition for the constant m (the eigenvalue of the projection of the momentum 
on the z axis) to be an integer appears due to the periodicity condition 

^(r, V5, z) = ^(r, + 2tt, z). 

Then for the radial wave function, ipmi'r) we then obtain the equation 

^9.(r9.V^„) + - ^^-^ - k-) i^m = 0. (99) 

We introduce the new radial coordinate 

P 



This is the transformation in (f7^ and (f7^ up to a constant. The radial equation is then 

-^dpipdpi^m) + - ^) = 0, (100) 

where 

2ME \m\ 

^ = - ^ ' = V ^ °- 

This is the Bessel equation. We solve it with the boundary condition 

^m|p=H/7 = 0, (101) 

which corresponds to the motion of an impurity inside the cylinder with an impenetrable 
boundary. The general solution of the Bessel equation ()100|) is 

i^ni = Cm-MPp) + dmKiPp), Cm, dm = COHSt, 

where J^, and N^, are Bessel and Neumann functions of order u [83J. The boundedness of 
the wave function on the axis of the cylinder requires d^, = 0- The constants of integration 
Cm are found from the normalization condition 

R 

drr\ipm? = 1- 
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Boundary condition ()10H) yields the equation for P 



U(5Rh) = 0, (102) 

which defines the energy levels of bounded states. It is well known that for real u > —1 
and -R/7, this equation has only real roots. Positive roots form an infinite countable set 
and all of them are simple [221 • This provides the inequality 

, 2ME , , 

We label the positive roots of Eqn ()102|) by the index n = 1,2,... (principle quantum 
number): j3 — > Pn{fn, a, 7, R)- Then the spectrum of bounded states is 

Er^ = ^ik' + Pl). (103) 
For large radii I3p^ 1 and (3p^ v we have the asymptotic form: 



npp V^^ 2 4 
As a result, we obtain explicit expression for the spectrum 



/3»=?f" + ^-T). (104) 



R \ 2a 4 

In the absence of defect, a = 1, 7 = 1, p = r, and the radial functions ipm are 
expressed through the Bessel functions of integer order u = \m\. In this case, the spectrum 
of bounded states depends only on the sizes of the cylinder. In the presence of the 
wedge dislocation, Bessel functions have a noninteger order in general. In this case, the 
spectrum of energy levels acquires dependence on the deficit angle 9 and the Poisson ratio 
a characterizing elastic properties of the medium. 

If the mass of impurity and vacancy is defined by integral pO|) . then M > for 
impurity and the energy eigenvalues are positive. For vacancy, M < and the energy 
eigenvalues are negative. In this case, the energy spectrum is not bounded from below, 
which causes serious problems for physical interpretation. It seems that one has to insert 
in the Schrodinger equation not the bare mass ()10|) but the effective mass, with the 
contribution of elastic stresses arising around a vacancy taken into account. This question 
presently remains unanswered. 

The presentation in the present section is close to that in j28j. In contrast to that 
paper, we use the elastic gauge for the metric. Therefore, our results depend not only on 
the deficit angle of a wedge dislocation but also on the elastic properties of the medium. 

Calculations of the energy levels of an impurity in the field of a wedge dislocation are 
actually equivalent to the calculations of bounded-state energies in the Aharonov-Bohm 
effect ^4j (see reviews jH3l|Hn]). The difference reduces only to changing the order of the 
Bessel functions, 

$ 



m 

$0 
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where $o = ^Tthc/e is the magnetic flux quantum and e is the electron charge. 

The considered example shows how the influence of the presence of defects is taken 
into account in the geometric approach in the first approximation. If calculations in 
some problem were performed in elastic media without defects, then to take the influence 
of defects into account we have to replace the flat Euclidean metric with a nontrivial 
metric describing the given distribution of defects. This problem may appear complicated 
mathematically because we have to solve the three-dimensional Einstein equations to find 
the metric. However, there are no principal difficulties: the effect of dislocations reduces 
to a change in the metric. 

15 Conclusion 

The geometric theory of defects describes defects in elastic media (dislocations) and defects 
in the spin structure (disclinations) from a single point of view. This model can be used 
for describing single defects as well as their continuous distribution. The geometric theory 
of defects is based on the Riemann-Cartan geometry. By definition, torsion and curvature 
tensors are equal to surface densities of Burgers and Frank vectors, respectively. 

Equations defining the static distribution of defects are covariant and have the same 
form as equations of gravity models with dynamical torsion. To choose a solution uniquely, 
one must fix the coordinate system. For this, the elastic gauge for the vielbein and the 
Lorentz gauge for the §0(3) connection are proposed. In the defect-free case the dis- 
placement vector field and the field of the spin structure can be introduced. Equations 
of equilibrium are then identically satisfied, and the gauge conditions reduce to the equa- 
tions of elasticity theory and the principal chiral SO(3)-field. In this way, the geometric 
theory of defects incorporates elasticity theory and the model of principal chiral field. 

In a certain sense, the elastic gauge represents the equations of the nonlinear elasticity 
theory. Nonlinearity is introduced in elasticity theory in two ways. First, the deformation 
tensor is defined through the induced metric 

instead of being defined by linear relation (0)). Then the stress tensor is given by an 
infinite series in the displacement vector. Second, Hook's law can be modified assuming 
nonlinear dependence of the stress tensor on the deformation tensor. Hence, the elastic 
gauge condition is the equations of the nonlinear elasticity theory where the deformation 
tensor is assumed to be defined through the induced metric and Hook's law is kept linear. 
A generalization to the nonlinear dependence of the deformation tensor on the stress 
tensor is obvious. 

As an example, we considered the wedge dislocation from the standpoint of the elas- 
ticity theory and the geometric theory of defects. We showed that the elasticity theory 
reproduces only the linear approximation of the geometric approach. In contrast to the 
induced metric obtained within the exact solution of the linear elasticity theory, the ex- 
pression for the metric obtained as the exact solution of the Einstein equations in the 
elastic gauge is simpler and is defined on the whole space and for all deficit angles. The 
obtained expression for the metric can be checked experimentally. 
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Two problems are considered as applications of the geometric theory of defects. The 
first is the scattering of phonons on a wedge dislocation. In the eikonal approximation, 
the problem is reduced to the analysis of extremals for the metric describing a given 
dislocation. Equations for extremals are integrated explicitly, and the scattering angle is 
found. The second of the considered problems is the construction of the wave functions 
and energy spectrum of impurity in the presence of a wedge dislocation. This requires 
solving the Schrodinger equation. This problem is mathematically equivalent to solving 
the Schrodinger equation for bound states in the Aharonov-Bohm effect |Hn]. The explicit 
dependence of the spectrum on the deficit angle and elastic properties of the medium is 
found in the presence of a wedge dislocation. 

The geometric theory of static distribution of defects can also be constructed for 
membranes, i.e., on the plane M^. For this, one has to consider the Euclidean version [SZ] 
of two-dimensional gravity with torsion [88-90] . This model is favored by its integrability 
[91-94]. 

The developed geometric construction in the theory of defects can be inverted, and 
we can consider the gravitational interaction of masses in the universe as the interaction 
of defects in elastic ether. Point masses and cosmic strings 1^ then correspond to 
point defects (vacancies and impurities) and wedge dislocations. In this interpretation of 
gravity, we have a question about the elastic gauge, which has direct physical meaning in 
the geometric theory of defects. If we take the standpoint of the theory of defects, then 
the elastic properties of ether correspond to some value of the Poisson ratio, which can 
be measured experimentally. 

It seems interesting and important for applications to include time in the considered 
static approach for describing motion of defects in the medium. Such a model is lacking 
at present. From the geometric standpoint this generalization can be easily performed, 
at least in principle. It is sufficient to change the Euclidean space M.^ to the Minkowski 
space M^'^ and to write a suitable Lagrangian quadratic in curvature and torsion, which 
corresponds to the true gravity model with torsion. One of the arising problems is the 
physical interpretation of the additional components of the vielbein and Lorentz con- 
nection that contain the time index. The physical meaning of the time component of 
the vielbein cq* OqU^ = is simple - this is the velocity of a point of the medium. 
This interpretation is natural from the physical standpoint because the motion of con- 
tinuously distributed dislocations means a flowing of the medium. In fact, the liquid can 
be imagined as the elastic media with a continuous distribution of moving dislocations. 
This means that the dynamical theory of defects based on the Riemann-Cartan geometry 
must include hydrodynamics. It is not clear at present how this can happen. Physical 
interpretation of the other components of the vielbein and the Lorentz connection with 
the time index also remains unclear. 

The author is sincerely grateful to I V Volovich for numerous discussions on the prob- 
lems considered in this paper. This work is supported by the RFBR (grant 05-01-00884) 
and the program of support for leading academic schools (grant AS-1542. 2003.1). 
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Appendix. Some differential geometry 



In the appendix, we briefly present the main facts from differential geometry and introduce 
the notation used in this review. The description is given in local coordinates, which is 
sufficient for our purposes. We recommend |i22J as a textbook on differential geometry. 

If the metric g and the affine connection F are given, then we say that the geometry 
is defined on a differentiable manifold M, dimM = m. We assume that all fields on the 
manifold are given by smooth C°°(M) functions except, possibly, some singular points, and 
do not stipulate that in what follows. We also assume that the manifold M is topologically 
trivial, i.e., diffeomorphic to the Euclidean space M™. 

In a local coordinate system x'^, /i = 1, . . . , m, the metric is given by a nondegenerate 
symmetric covariant second rank tensor g^y[x), which defines the scalar product of vector 
fields X = X^d^, Y = Vd^ 

(X, Y) = X^Y'^g,,, g^, = g,,, detg,, ^ 0. (105) 

In general, the scalar product may be not positive definite. If the scalar product is positive 
definite, we say that a Riemannian metric is given on a manifold. By definition, the metric 
is a covariant tensor field, i.e., it transforms under coordinate changes — > x^ (x) by the 
tensor law 

dx'^ dx" 

This means that the scalar product of vector fields (X, Y) is a scalar field. In a similar 
way, contractions with the metric g^^ and its inverse g^^ , g^'^gup = ^p-, allows one to build 
scalar fields from higher-rank tensors or to lower their rank. The metric is also used for 
lowering and raising of tensor indices. 

An affine connection on a manifold in a local coordinate system is given by the set of 
coefficients T ^^^^x) that transform under the diffeomorphisms as 

^ pi dx^ dx"^ ^ pdxP' ^ d'^xP dx^' (106) 
dx^' dx'^' dxP dxf^' dx'^' dxP 

These coefficients do not constitute a tensor field because of the presence of an inhomoge- 
neous term in (jlOfij) . An affine connection on a manifold defines covariant derivatives of 
tensor fields. In particular, covariant derivatives of a vector field and 1-form A = dx^A^ 
have the form 

V^X^ = 9^X'^ + XT^/, (107) 
V^A, = d^A, - T^/Ap. (108) 

The covariant derivative of a scalar field coincides with the partial derivative V/i(/J = dp^ip. 
The covariant derivative of higher-rank tensors is built in a similar way and contains one 
term with the plus and minus sign for each contravariant and covariant index, respectively. 
One can easily check that the covariant derivative of a tensor of an arbitrary type (r, s) 
is a tensor field of type (r, s + 1), i.e., it has one additional covariant index. We note that 
the covariant derivative of a product of tensor fields may contain arbitrary contractions 
of indices. For example, 

dpiX^A,) = VpiX-'A,) = {V.xnA, + X'^iV^A,). 
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Because the inhomogeneous term in ()106|) is symmetric in indices /x' and v' the anti- 
symmetric part of the affine connection 2T\p^yf forms a tensor field of type (1,2), which 
is called the torsion tensor 

V = V-r,/. (109) 

In general, the connection T has no symmetry in its indices and is not related to 
the metric g^y in any way because these notions define different geometric operations on 
a manifold M. We then say that the affine geometry is given on M. We emphasize that 
the metric and the affine connection are defined arbitrarily and are completely indepen- 
dent geometric notions. Therefore, in the construction of physical models, they can be 
considered independent fields having different physical interpretations. 

If the affine geometry is given on a manifold, then we can construct the nonmetricity 
tensor Qfj_yp that is by definition equal to the covariant derivative of the metric: 

By construction, the nonmetricity tensor is symmetric with respect to the permutation 
of the last two indices: Qfj,up = Qp,pu- We note that we need both objects, the metric and 
the connection to define the nonmetricity. 

Thus, for a given metric and connection we constructed two tensor fields: the torsion 
and the nonmetricity tensor. We prove that for a given metric, torsion, and nonmetricity 
tensor we can uniquely reconstruct the corresponding affine connection. Equation (jllOj) 
can always be solved for the connection F. Indeed, the linear combination 

y^gup + ^uQpp - ^pQ^iv 
yields the expression for the affine connection with all lowered indices: 

^ pup ^ pu Qap 2^^^^'^'^ ~^ (^uQpp dpQpy) + — (Tp_i,p Tupp + Tppy) 

~i~ i2^Q^iup ~\~ Qupp Qppu)- (H-'-) 

The right hand side of this equality is symmetric in the indices /i and v except one term 
Tpypjl^ and this is in accord with the definition of the torsion tensor ()109|) . Thus, to define 
the affine geometry on a manifold M, it is necessary and sufficient to define three tensor 
fields: metric, torsion, and nonmetricity. We stress once again that all three object can be 
specified in a completely arbitrary way, and they can be considered different dynamical 
variables in models of mathematical physics. 

It is easy to compute the number of independent components of connection, torsion, 
and nonmetricity tensors: 

71-^3 PI vn?ijn-\) ^ m2(m + l) 



pvp\ 



This shows that the total number of independent components of the torsion and non- 
metricity tensors equals the number of components of the affine connection. 
We now consider particular cases of the affine geometry. 
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In the attempt to unite gravity and electromagnetism, H. Weyl considered the non- 
metricity tensor of a special type [HH] 

Qm-p = W^g^p, (112) 

where is the Weyl form identified with the electromagnetic potential (here, torsion 
was assumed to be identically equal to zero). We say that the Riemann-Cartan-Weyl 
geometry is defined on a manifold if there are given a metric, torsion, and nonmetricity 
of special type (|112|) . 

If the nonmetricity tensor is identically zero, Q/j^up = 0, but the metric and torsion are 
nontrivial, then we say that the Riemann-Cartan geometry is given on a manifold. As a 
consequence of fll 1 1|) . the affine connection is defined uniquely by the metric and torsion 
in this case: 

1 1 
Tpyp = —{dp^g^jp + d^g^p — dpgp^) + —{Tp^p — T^pp + Tpp^). (113) 

This connection is called metrical because the covariant derivative of the metric is iden- 
tically equal to zero: 

"^pQup = dpg^p - Tp^^g^p - Tpp^g^^ = 0. (114) 

The metricity condition provides commutativity of covariant differentiation and the raising 
and lowering of indices. 

If the torsion tensor identically vanishes, Tpi^f = 0, and nonmetricity have special form 
()112j) . then we say that the Riemann-Weyl geometry is given. 

If both nonmetricity and torsion tensors are identically equal to zero, Qpup = 0, 
Tpup = 0, and the metric is nontrivial, then we say that Riemannian geometry is given on 
a manifold. In this case, the metrical connection is symmetric with respect to the first 
two indices and uniquely defined by the metric 

f A^i^p = ^{dpgup + dugpp - dpgpu). (115) 

This connection is called the Levy-Civita connection or Christoffel symbols. 

We use the tilde to denote the geometrical objects constructed only for the metric 
and zero torsion and nonmetricity, i.e., in a (pseudo-) Riemannian geometry. The prefix 
pseudo- is used when the metric is not positive definite. 

From the expression for Christoffel symbols ()115|) . we see that 

dpgup = ^fiup + ^f,pu- (116) 

This means that for the Christoffel symbols to be zero in some coordinate system it is 
necessary and sufficient that the metric components be constant in these coordinates. In 
another coordinate system, they may be nontrivial because Christoffel symbols are not 
components of a tensor. For example, the Christoffel symbols for Euclidean space are 
zero in the Cartesian coordinates but differ from zero, e.g., in spherical or cylindrical 
coordinate systems. 
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In the case where the nonmetricity tensor and torsion are identically zero, and there is 
a coordinate system in which the metric is equal to the diagonal unit matrix in the neigh- 
borhood of every point, and, consequently, the Christoffel symbols are zero, the geometry 
is called locally Euclidean. The corresponding coordinate system is called Cartesian. 

We say that a vector field X = X^d^^ is transported parallel along a curve x^(t), t G R, 

if 

x^'V^X^' = 0, 

where x'^ = dx^ /dt is the tangent vector to the curve. Multiplying this equation by dt we 
obtain 

= dx^d^X^" = -bx^'Y^p^XP, 

where 6x^ = dtx^ . We say that under the parallel transport from a point x'^ to a neigh- 
boring point x'^ + 5x'^ , the vector X^" acquires a differential 5X^, which is linear in 6x" 
and the components of the vector field. In a similar way, we define parallel transport 
of arbitrary tensor fields along a curve The result of parallel transport between 

two points depends in general on a curve connecting these points. We note that parallel 
transport along a curve is defined only by an afiine connection and has no relation to a 
metric. 

There are two types of distinguished curves x^{t) in an affine geometry: geodesies and 
extremals. A geodesic line is a curve such that a vector tangent to it remains tangent 
under parallel transport along the curve. With the parameter t along the curve chosen 
canonically, a geodesic line is defined by the system of ordinary nonlinear differential 
equations 

x^' = -T^/x'^xP. (117) 

Although a geodesic line is defined only by the symmetric part of the afiine connection 
T^pujp, the latter depends nontrivially on torsion and nonmetricity. Explicit expression 
(|TTT|1 yields 

The second type of distinguished curves in an affine geometry are extremals or lines 
of extremal length connecting two points. These lines are exclusively defined by the 
metric. With the parameter along the curve chosen canonically, extremals are defined by 
an equation similar to pi7|) . 

x^ = -f^/x^x". (118) 

However, we now have not a general affine connection on the right-hand side but Christof- 
fel symbols constructed only from the metric. We see that geodesies and extremals are 
in general different curves on a manifold. In a Riemann-Cartan geometry, geodesies and 
extremals coincide if and only if the torsion tensor is antisymmetric in all three indices. 
In Riemannian geometry, geodesies and extremals always coincide. 

The primary role in differential geometry is played by the curvature tensor of the affine 
connection, which arises in different contexts. In local coordinates, it is defined by the 
expression 

R^.up'' = dpT,/ - r^p^r.A" - (/i ^ ly), (119) 
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where the parenthesis (/x ^ u) denote the previous terms with the permuted indices fj, 
and u. It can be easily verified that curvature ()119|) is indeed a tensor field. We note that 
the curvature tensor has no relation to a metric and is defined entirely by the connection. 
In an affine geometry, the curvature tensor with all of its indices lowered has no symmetry 
under permutations of indices except antisymmetry in the first two indices. 

Contraction of the curvature tensor in the two indices yields the Ricci tensor R^^, = 
Rfipv^ 1 which is also defined entirely by the connection. In Riemannian geometry, the 
Ricci tensor for the Levi-Civita connection is symmetric with respect to the permutation 
of indices: R^j_i, = Ru^- In a Riemann-Cartan geometry, this symmetry is generally absent 
for a nonzero torsion tensor. 

If a metric is also given on a manifold, then we can construct the scalar curvature 
R = R^^g^^r 

The affine connection is called locally trivial if in the neighborhood of every point one 
can choose the coordinate system in which all components of the connection are zero. 
Now we formulate two important theorems. 

Theorem 1. For the local triviality of the affine connection it is necessary and sufficient 
that its torsion and curvature tensors are equal to zero on M. 

The proof of this theorem is reduced to an analysis of the transformation rule for the 
affine connection components in ()106|1 . If in the new coordinate system the symmetric part 
of the connection components is equal to zero, then the system of differential equations 
for the transition functions appears. The local integrability condition for this system of 
equations is provided by the equality of the curvature tensor to zero. 

A more thorough statement is proven in [99j. 

Theorem 2. // for a given affine geometry the torsion, nonmetricity, and curvature 
tensors are equal to zero on the whole manifold, then this manifold is isometric either 
to the whole (pseudo-) Euclidean space or to the quotient space /Gs, where G is a 
discreet transformation group acting freely. 

In the last theorem, we suppose that both the metric and affine connection are given 
on a manifold. 

The Riemann-Cartan geometry defined by the metric and torsion provides the basis 
for the geometric theory of defects. In the analysis of such models, it is more convenient to 
introduce the Cartan variables: the vielbein and SO(m) connection instead of the metric 
and torsion. We assume here that a metric on a manifold is positive definite. If a metric 
were not positively definite then the §0(]9, q) connection would appear, where p + q = m. 

For a Riemannian metric, the orthonormal vielbein e^*(x), i = 1, . . . , m, is defined by 
the system of quadratic equations 

Qfiu = e^^'e^Sij, (120) 

where Sij is the Kronecker symbol. This relation uniquely defines metric for a given 
vielbein. Conversely, system of equations (|12Up defines the vielbein for a given metric 
up to local SO(m) rotations. We distinguish Greek and Latin indices because different 
transformation groups act on them. From definition ()120|) . we see that dete^* 7^ 0. 
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Components of the inverse vielbein e^j, e^ie^ = 6^, can be considered components of m 
orthonormal vector fields Cj = e'^j(9^ on a manifold M with respect to the metric g 

Components of the vielbein e^* define m orthonormal 1-forms e* = dx^e^^ on a manifold. 

It is known that any manifold M can be equipped with a Riemannian metric (see, for 
e.g., |lUUj ). At the same time, the global existence of a vielbein provides, in particular, 
an orientation of the manifold M. Thus, a vielbein may exist globally only on orientable 
manifolds. There are also other topological restrictions for the global existence of a 
vielbein that we do not discuss here. 

The coordinate basis of the tangent space in every point of M is called holonomic. 
Components of tensors of an arbitrary type may also be considered with respect to the 
unholonomic basis Cj of tangent spaces defined by the vielbein. For example, a vector 
field has the components 

X = X'^d^ = X'a, X' = x^e;. 

We always assume that the transformation of Greek indices into Latin ones and vice versa 
is performed with the help of the vielbein. 

We now define the SO(m) connection uj^i^{x). This is the name of the connection of 
the principal fiber bundle with the structure Lie group §0(m) and the base M. If the 
Riemann-Cartan geometry is given on M, and the vielbein is defined, then we can define 
the SO(m) connection by the relation 

V^.eJ = d^ej - T^/e/ + e,^u^/ = 0. (121) 

We see that under coordinate changes, the components of tu^/ transform as a covector 
field in the Greek index. Thus, they define a 1-form on M. For a given vielbein, this 
relation provides a one-to-one correspondence between components of the connections 
Tfj^iyP and uj^j^. The §0(m) connection defines covariant derivatives for components of 
tensor fields relative to an unholonomic basis. For example, 

V^X' = d^X' + X^LU^/, V^Xi = d^Xi - iu^^X,. (122) 

The covariant derivative is naturally defined for a tensor field having both Greek and 
Latin indices. Taking the covariant derivative of Eqn p20|) leads to the antisymmetry of 
the components cj^*-' = —uj^^^. This means that the 1-form dx^u^j"^ takes values in the 
Lie algebra so(m), and this indeed corresponds to the SO(m) connection. 

In the general case of an affine geometry, the Cartan variables can also be defined by 
relations ()12U|) and p2H) . For this, we have to replace the Kronecker symbol on the right- 
hand side of ()120j) with an arbitrary nondegenerate symmetric matrix rjij. In this case, 
relation (|12H) defines a linear GL(m, R) connection. In the Riemann-Cartan geometry 
with the Lorentzian-signature metric, we would have the Lorentz §0(1, m — 1) connection. 

We consider local rotations with a matrix S/{x) G SO(m) to show that the components 
a;^/ do define an SO(m) connection in the Riemann-Cartan geometry. By definition, the 
components of vector fields and 1-forms transform under local rotations according to the 
rule 

X'' = X^Sj\ X'. = S-\^Xy (123) 
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In order that covariant derivatives ()122|) have the tensor transformation rule under local 
rotations, it is necessary and sufficient that the components of the SO(m) connection 
transform according to the rule 

u'J = S-\^u,i}Si' + d^S-\''Sk'. (124) 

This is the transformation law for the SO(m) connection, indeed. The same transfor- 
mation law follows from definition ()12H) if the vielbein e^* transforms as a vector with 
respect to the index z, and the Christoffel symbols remain unchanged. 

Of course, we can introduce the metrical and SO(m) connections on M in an inde- 
pendent way. If we require afterwards that the SO(m) connection acts on tensor field 
components relative to an unholonomic basis defined by the vielbein, we obtain a one-to- 
one correspondence between the components of the connections ()121|) . 

We now express the components of the metric connection T ^j^p through the vielbein 
e^' and §0(m) connection tu^j* with the help of relation ^V21\ji and substitute them into 
the definitions of torsion (jl(J9p and curvature (|119|) . As a result, we obtain expressions 
for torsion and curvature in terms of the Cartan variables 

T^u = dpe.' - - ^ (125) 

Rpvj = d^^ujyj' - ujpj'^ujyk - (Ai ^ v), (126) 

where 

Torsion p25|) and curvature ()126|) are 2-forms on the manifold M with the values in the 
vector space and Lie algebra so(m), respectively. 

An SO(m) connection uo^j^ is called locally trivial if every point has a neighborhood 
containing this point such that it has the form 

= d^S-^,^Sk\ (127) 

Obviously, after local rotation with the matrix 5*"^/, all components of the connection 
become zero. This connection is also called a pure gauge. One can easily verify that the 
curvature of a locally trivial connection identically vanishes. The inverse statement is also 
valid. 

Theorem 3. An §0(m) connection is locally trivial if and only if its curvature form is 
identically zero on M. 

This theorem is valid for any structure Lie group. The proof is reduced to an analysis 
of transformation rule ()124|) . If the left-hand side of this equation is zero, then we have 
a system of equations for the field S/{x) for which vanishing of curvature tensor is the 
local integrability condition. 

The space with zero curvature tensor Rpuj^ = is called the space of absolute paral- 
lelism because the parallel displacement of a vector does not depend on the path connect- 
ing two fixed points of the manifold. 

We can perform a local SO(m) rotation for the locally trivial §0(m) connection such 
that it becomes zero: cj^/ = 0. Then, the zero-torsion equality becomes 

dpej - d^e.' = 0. 
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This equation is the local integrability condition for the system of equations 

d^^y' = e/ (128) 

for m functions y'' {x) . A solution of this system of equations yields the transition functions 
to Cartesian coordinates. Thus, the equalities of curvature and torsion tensors to zero 
are the necessary and sufficient conditions for the existence of the fields S/{x) and y''{x), 
i.e., the existence of such a local rotation and coordinate system where the connection 
vanishes and the metric becomes Euclidean. We note that equality of torsion tensor to 
zero alone is not enough for the existence of a Cartesian coordinate system. 

Three-dimensional space is considered in the geometric theory of defects. We make 
two remarks concerning this. In lower dimensions, the algebraic structure of the curva- 
ture tensor with all lowered indexes becomes much simpler. In two dimensions the full 
curvature tensor in the Riemann-Cartan geometry is in the one-to-one correspondence 
with its scalar curvature: ^ 

Rijkl = -^{SikSjl — SiiSjk)R. 

In three dimensional space, the full curvature tensor is in the one-to-one correspondence 
with its Ricci tensor 

Rijkl — ^ikRjl — ^ilRjk — ^jkRil + ^jlRik — ^{^ik^jl — 5il5jk)R- (129) 

These formulas are also correct for nonzero torsion. 

Concluding this appendix, we write the identity that is valid in the Riemann-Cartan 
geometry in an arbitrary number of dimensions: 

i?(e, uj) + ^TijkT''^ - ^TijkT"'^ - T/T' - -d^{eT'') = ^(e), e = det (130) 

The Riemannian scalar curvature on the right-hand side of this equality is constructed 
only from the vielbein for zero torsion. 
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